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SECTION 1.0 
INTRODUCTION 

1.1 Background 

Wind shear j particularly at the lower altitudes in the terminal area, 
has been identified as being hazardous to aircraft operations. Accurate and 
reliable wind profiles are required for use in fast time and manned flight 
simulation studies aimed at fully defining and understanding the wind shear 
hazard. This report describes wind speed profiles developed for the above 
simulation to improve the safety and reliability of operations in the 
terminal area. A comprehensive set of wind profiles and associated wind 
shear characteristics which encompass many of the wind shear environments 
potentially encounterable by aircraft in the terminal area have been modelled, 
subject to the data available. For the purpose of this effort, wind shear is 
defined as significant changes in wind speed and/or direction up to 500 m 
above the ground which may adversely affect the approach, landing, or takeoff 
of an aircraft. The wind shear is mathematically modelled and the mathe^ 
matical scenarios (environment) are presented in a format for direct applica- 
tion to wind shear hazard/flight simulation studies. 

1.2 General Description of Wind Shear Models 

A survey of existing wind shear data and mathematical models which are 
comprehensive, of sufficient spatial extent, and detailed enough for the 
development of mathematical models of wind shear is reported in Reference 
[1-1]. This reference reviews the state of the art and describes the scale 
and duration of frontal wind shear, neutral and stable boundary layer wind 
shear, and thunderstorm wind shear. The physics of these atmospheric phenom- 
ena with which significant wind shear is associated is also outlined in this 
reference and is not described in the present report. Selected data sets 
reported in the literature and described in Reference [1-1] have been chosen 
for the three wind shear conditions considered, i.e., 1) thunderstorm wind 
shear, 2) neutral and stable boundary layer wind shear and 3) frontal wind 
shear. In all three cases, the most comprehensive set of wind shear data 
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reported in the literature have been discretized to a grid system and a 
computer program lookup with interpolation capabilities developed. The 
spatial variation of wind speed, both horizontally and vertically, is obtained 
from the programs when called by a main simulation program. The location of 
data measurement, how the data were obtained, the time of observation and 
the extent of the data up to 500 m are outlined for each of the wind shear 
conditions considered. From these data sets, a description of steady state 
wind which, in the thunderstorm and frontal cases, is effectively a quasi- 
steady state wind having a TO sec averaging period are described. For the 
stable and neutral boundary layer case, the wind speeds are mean values 
which are defined as a 10 min average or greater. For the neutral and 
stable boundary layer, spatial dependence of the wind only in the vertical 
direction, z, is considered and the wind field in these cases is assumed 
homogeneous in horizontal extent. For the thunderstorm and frontal wind 
fields, both vertical, z, and horizontal, x, spatial dependence is considered. 
The data available, however, represents wind fields in a vertical plane 
sliced through the storm and spatial variation in the direction lateral to 
the motion of the storm is not considered. All three velocity components, 
that is, longitudinal velocity, W , lateral velocity, W , and vertical 
velocity, W_, are defined at each point in the vertical spatial plane, 
however. 

Since the high frequency component of the variation in the wind has 
been averaged out of the data sets used to define the mathematical form of 
the wind shear models, turbulence simulation techniques are provided from 
which a turbulent fluctuation can be generated and superimposed upon the 
quasi-steady state wind fields. For the thunderstorm and frontal models, 
the frequency content of the turbulent simulation must contain all fluctua- 
tions greater than 0.1 Hz. For the neutral and stable boundary layer the 
turbulence model must contain fluctuations of all frequencies. 

Statistical models which will allow estimates of the probability of a 
given shear magnitude and the frequency of occurrence are developed to the 
extent current data permits. A statistical description of the wind fields 
are required to establish meaningful magnitudes of wind shear to assure that 
the values used in flight hazard studies are realistic and might actually be 
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encountered in the real world. Statistical data, however, is very limited 
in this regard and it was necessary in most cases to establish only crude 
estimates of the probability of wind shear and the risk of exceeding a 
prescribed value. 

1.3 Organization 

The mathematical models of wind speed developed for thunderstorms, 
stable and neutral boundary layers, and the synoptic frontal storms are 
reported individually in Sections 2, 3 and 4, respectively. Each section 
describes the nature of the data set used to develop the table lookup programs 
and the theoretical turbulent simulation models recoranended for use with 
the quasi-steady winds. Statistical models are also described and estimates 
of probabilities and risk of exceedance given. In the appendix, tables of 
the wind fields for easy engineering application are given. Also, numerous 
illustrative plots of typical wind speed profiles encountered along conven- 
tional flight paths are provided for visualization of the type and magnitude 
of wind shear potentially encounterable in the atmospheric conditions con- 
sidered. Finally, a careful listing of the computer code and the associated 
data set are described. The computer programs are givenM'n the form of 
subroutines which can be irwnediately coupled with a fast time computer 
program for calculating flight through wind shear hazards or for direct 
application to flight simulation studies. This material has been placed in 
the appendices because of its bulkiness and because, although it directly 
pertains to the engineering application of this report, it is not immediately 
relevant to the text of each section. However, the appendices have been 
numbered corresponding to the section of the report to which they refer. 
Appendix 1 contains the nomenclature while Appendices 2, 3 and 4 contain the 
tables, illustrations and computer code for the respectively numbered 
sections. Finally, each section contains recommendations and guidelines for 
engineering applications. The range of applicability, engineering interpre- 
tation, and examples of how to use the data are given. 

The report concludes with a surrenary section which summarizes the results 
of the effort, presents conclusions, and provides reconmendations concerning 
future efforts. A discussion of data deficiencies and the associated weak- 
nesses and strengths of the model are also described in the summary section. 
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SECTION 2.0 


THUNDERSTORM WIND SHEAR 

A model of wind shear associated with thunderstorms has been developed 
from the data of Goff [2-1], These data consist of measurements of longi- 
tudinal, lateral and vertical wind speeds at various levels of a 500 m 

tower. Although the data are measured along a vertical line in the atmo- 
sphere, they are projected horizontally to form a vertical plane using the 
concept of frozen turbulence or Taylor's hypothesis. Thus, the wind shear 
phenomena associated with thunderstorm models can be described as a two- 
dimension, spatial wind field. 

The data of Goff [2-1], after careful review (see Reference [2-2]), 

were selected as being the best data available to construct a quantitative 

model of thunderstorm wind shear that provides both vertical and horizontal 
wind shear values. As noted, the model is restricted to two-dimensional 
space. The extensive survey reported in Reference [2-2] indicates that 
there are no three-dimensional data sets available nor any theoretical 
models associated with thunderstorms which would allow a three-dimensional 
simulation to be carried out. Therefore, it was necessary to restrict the 
simulation of thunderstorm wind shear to two-dimensional wind fields. The 
data, however, do include all three velocity components of the wind vector 
in the plane swept out by the 500 m tower as the storm passes. 

2.1 Data Source 

Thunderstorm wind shear data are presented in the form of longitudinal 

(W ), lateral (W ) and vertical (W ) wind speed components in a vertical 
X y z 

plane for 20 thunderstorm cases. Data for these 20 storms were measured 

during the months of May through June over the period of 1971 through 1973 

with the WKY-TV/NSSL 481 m meteorological tower, Norman, Oklahoma. Time 

histories of the wind speed are converted to horizontal spatial distributions 

using Taylor's hypothesis (i.e., x = Wt). Ten second averaged values of 

wind speed components are provided in the form of isotach maps for W , W 

X y 

and W^, respectively. These data were converted to a 41 x 11 point grid 
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system as illustrated in Figure 2-1. The data for the 20 thunderstorm cases 
having been discretized on the grid system were stored on magnetic tape and 
a computer lookup routine developed for interpolating the W^, and wind 
speeds for any location within the x-z plane, tabulated data for W , W 
and On the 41 x 11 point grid system are given in Appendix 2A; a pictorial 
description of the wind speed profiles and streamline patterns are given in 
Appendix 2B, and a computer program which stores the data on disks and 
carries out the table lookup with the option of superimposing turbulencO is 
described in Appendix 2C. 

The 20 thunderstorm cases for the purpose of this report are assigned 
numbers 1 through 20 which correspond to Goff's identification symbols A 
through T. All tabulated and illustrated data in the appendix contain both 
the numerical and the alphabetical identification. 

2 . 2 Data Process i nq 

For each thunderstorm, a 10 min record was taken and the data were 
averaged over 10 sec intervals throughout the 10 mins. These data were then 
fitted to a regular array of a 60 x 10 grid by Goff. Although the actual 
measurements were not uniformly spaced, the data as presented for public 
distribution were uniformly spaced on the 60 x 10 grid. Each 10 min record 
consists of 60 10 sec averaged data sets representing a cross section through 
the storm. Streamline patterns as well as isotach maps of W , W and W 

X Jr ^ 

were also drawn from the data and are given in Reference [2-1]. 

The data presented in Reference [2-1] were evenly spaced over 10 50 m 
intervals in the z-direction. Values assigned to each interval were deter- 
mined by linear interpolation from seven levels of measurements for the 
1971, 1972 data and from only three levels of measurement for the 1973 and 
1974 data. For the x-direction, 60 intervals were utilized where a Galilean 
transformation from time to space was utilized. The horizontal extent of 
each cross-section is therefore a variable and is equal to the frontal speed 
of the thunderstorm times the 10 min averaging period. 

The coordinate system utilized in Reference [2-1] and also in this 
report defines the x-direction as being measured positively in the direction 
of frontal motion and the y-di recti on as being measured perpendicular to the 
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frontal axis utilizing the right-hand rule, i.e., positive is measured 
into the plane of the paper. The W components of the wind speed are posi- 

A 

tive in the direction of frontal motion, i.e., positive from left to right on 
the contour maps and the vertical wind speed, W^, is positive in the upward 
direction. 

The data presented by Goff [2-1] are the longitudinal wind speed, W^, 

A 

relative to the motion of the storm. In this report, however, the data are 
generally presented relative to the fixed frame of reference attached to the 
ground. 

In view of the 10 sec averaging period utilized to reduce the data, all 
fluctuations in wind speed of frequency higher than 0.1 Hz have been filtered 
out of the data. Thus, the data contain variations in wind speed having fre- 
quencies of 0.1 Hz or less. These low frequencies are expected to be those 
to which aircraft motion is most sensitive. However, to assure correct 
simulation where high frequency wind components are significant, it is 
recomnended that a turbulent fluctuation be superimposed upon the quasi - 
steady wind fields. Section 2.4 discusses and recommends a turbulence model 
for use with the thunderstorm gust front data. 

The following section. Section 2.3, describes the quasi-steady wind 
speed profile grid system and presents working data for wind shear hazard/ 
flight simulation studies. The wind speed is referred to as quasi-steady in 
view of the fact that it is averaged over 10 sec time intervals and thus 
contains departures from the mean wind speed which is generally averaged over 
a 10 min period or greater. 

2.3 Quasi-Steady Wind Speed Profile Grid System 

The wind fields of Goff [2-1] were fit to a 41 x 11 point grid system. 
The grid system as illustrated in Figure 2-1 was superimposed on the contours 
and wind speeds at each grid point tabulated and punehed onto computer cards. 
The data were later stored on magnetic tape. 

The grid system is numbered from the left-hand bottom corner. The 
numbers increase from left to right in the positive x-di recti on of the 
original data and from bottom to top in the positive z-direction of the 
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original data. The wind speeds are given in units of meters per sec, 
m s“^, with W being positive in the direction of frontal motion, W being 

J\. £m 

positive upward and W being positive into the plane of the paper. 

The wind speed W is stored on tape as the wind speed relative to the 
storm motion. The stored values of wind speed are therefore the values an 
observer would measure moving along with the storm. To convert to the earth 
frame of reference the mean motion of the storm must be added to the longi- 
tudinal wind speed. The next section describes tables and graphical illustra- 
tions of the wind speed profiles. In these cases the mean motion has been 
added . 

2.3.1 Tables of Wind Speed 

Tables of thunderstorm wind speed are given in Appendix 2A. The upper 
portion of the table covers grid stations 1 through 21 and the lower table 
covers grid stations 21 through 41. (Note column 21 is repeated for symmetry 
and clarity of presentation.) 

The thunderstorm case numbers designation for this report are listed at 
the top of each table. The letters in parentheses and the following series 
number corresponds to the thunderstorm designation given by Goff [2-1]. 

Also listed at the top of the table is the frontal speed, W^, and the 
horizontal length scales for the given wind record. The horizontal extent of 
each data set varies because of the data reduction procedure. Hence, the 
length of field, L, in kilometers and the horizontal grid spacing. Ax, are 
specified on each table. The vertical extent of each field is taken as 500 m 
with 50 m vertical grid spacing. 

2.3.2 Illustrations 

Illustrations of the longitudinal, lateral, and vertical wind speeds 
encountered along a 3® glide slope through each thunderstorm are provided in 
Appendix 2B. The glide slope is adjusted to terminate at the lower left-hand 
corner of each streamline plot as illustrated. The ordinate in the wind 
speed profiles is height, z, nond i mens ionali zed with the length, H = L tan 3°, 
where L is the length of the wind field. Each profile is the wind seen by an 
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airplane traveling along the flight path line drawn across the streamlines as 
shown in the upper figures. The flight paths are drawn to terminate in the 
lower left-hand corner always. Note that the streamlines plotted are 
relative to the speed of the front which for reference purposes is indicated 
with the vertical dashed line on the horizontal wind speed profile. The wind 
speed profiles are relative to the fixed earth frame of reference. 

2.3.3 Computer Program for Computing Thunderstorm Mind Speed Profiles 

A computer program has been developed which for given input x, z 
computes the longitudinal, lateral, and vertical wind speeds at that position. 
The six velocity gradients 3W /3x, 9W /9z, 9W /9x, 9W /9z, 9W /9x, and 9W /9z 

X X ^ ^ jr Jf 

are also output by the program. A detailed description and user's instruc- 
tions for the computer program are provided in Appendix 2C. The computer 
program also has the option of calling for turbulence which is super- 
positioned upon the thunderstorm wind field. A non-Gaussian turbulence model 
based on the technique of Reeves, et al . [2-3] is employed for turbulence 
simulation. 

It is intended that this computer program can be used as a subroutine 
for direct application to fast time flight studies or for fast time computer 
simulation of aircraft flight through thunderstorms. As an example of the 
application of the computer program, the reader is referred to Reference [2-4] 
in which this computer program has been used to study aircraft dynamics in 
thunderstorms and Reference [2-5] where the data were used for manned flight 
simulator programming. 

2.4 Turbulence in Thunderstorms 

Measurements of the power spectral density function for turbulence in 
thunderstorms has been reported as early as 1962 by Steiner and Rhyne [2-6]. 
Their data was measured only over the approximate reduced spatial frequency 
range of 0.004 to 0.4 rad ra”^ . The theoretical von Karman spectrum follows 
the data in this frequency range very well as demonstrated by Houbolt, 
et al . [2-7], in Figure 2-2a. The Dryden spectrum, on the other hand, does 
not compare as well with the data. However, with the properly chosen length 
scale a reasonable fit is obtained (see Figure 2-2b). All reported data. 
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(a) von Karman Model 



Reduced Frequency, Radians/m 
(b) Dryden Model 


Figure 2-2 Measured and Fitted Spectra for Thunderstorm . L'vlntegral 
Length Scale [2-7] 
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however, were measured in the altitude range of 12 to 8 km and moreover they 
do not extend to low enough frequencies to illustrate at which point the knee 
of the turbulence spectrum curve occurs. The knee of the power spectral 
density function is important in many design applications (see Reference 
[2-8]) and particularly to the empirical value of length scale, L, used in 
the analytical models of the power spectral density function. 

In Reference [2-7] a comparison of the power spectral density function 
of severe storms with that of cumulus clouds and clear air turbulence is 
given (see Figure 2-3). One can see from Figure 2-3 that the turbulent 
spectra for severe storms behaves very similar to that of cumulus clouds and 
clear air turbulence with the only major difference being higher values of 
<|)(fi) which indicates higher turbulence intensity. Again, these data are 
measured at very high altitudes and probably do not include effects due to 
the presence of the ground. 

Reference [2-7] computes two values of turbulence intensity. One is 
computed directly from the experimental data by integrating the power 
spectral density function over only the limits of the frequency range in 
which the data have been measured, i.e.. 




(})(fi)dfi 




(2-1) 


where u and i designate upper and lower, respectively. Spatial frequency 
designated by is related to cyclic frequency used in Section 3-4 by 
n = Wfi/2ir. An alternate method of measuring the turbulence intensity is to 
compute its value from the autocorrelation function evaluated at zero. 


= R(o) 


oo 

r 


(J)(Ji)dQ 

o 


( 2 - 2 ) 


It is found that based on the Dryden model the two measures of turbulence 
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Figure 2-3 Typical Power Spectra of Vertical Component of Turbulence 

Measured in Clear Air, Cumulus Cloud, and Thunderstorm [2-7] 



intensity are related by 




(2-3) 


where and represent the upper and lower limits, respectively, of the 
range of the turbulent spectrum measurements. Also appearing in this equa- 
tion is the length scale, L, which is an extremely important parameter in 
fitting the analytical expressions for the turbulence spectra to the data 
(see Reference [2-9]). 

Although the von Karman spectra tends to fit the data better as 

illustrated in Figure 2-2, the Dryden spectrum is conventionally used because 

! 

of its rational form. 


‘■w^'w 

X X 

f \ 

1 

7T 

1 + (L. , 
w.. 


[1 f (L !1)27 


(2-4) 



2tt 


[1 + 3(L n)^] 

z 

[1 + (L S1)2]2 

z 


Evaluation of the Dryden spectrum in a turbulence simulation routine 
requires knowledge of the length scale, and of the turbulence intensity, 

3. Houbolt, et al. [2-7] recommend values of L = 1036 m and gives ai values 
for the vertical fluctuations in the range of 4.88 m s”^ to 2.27 m s"^ and 
for the lateral fluctuations of 4.69 m s”^ to 2.69 m s”\ Longitudinal 
values of o were not measured. Adjusting oi to a from Equation 2-3 gives 
values of approximately 10.20 > > 4.75 m s"^ and 9.82 > > 5.63 m s V. 
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These values cannot be used directly, however, because the ground is 
expected to have a strong effect on the turbulence length scale and intensity. 


No actual data for L and o nor for how they vary with height in a 
thunderstorm below 500 m appears to be available in the literature. Barr, 
et al. [2-10] postulate a decrease in length scale and intensity at low 
altitudes. Preliminary results from the analytical model of Lewellen, 
et al. [2-11] predict increasing values of and Oyi- near the ground (see 
Figure 2-4). In Section 3.4 it is recommended that near the ground the ratio 
of and be determined from Figure 3-5. The relationships are 

/a =[0.177 + 2.74 X 10“^ z]"°*^ 

X z 


and 


a /a 

w/ 


[0.583 + 1.39 X 10"^ Z']"°*® 


(2-5) 


and for lei,ck of a better model this relationship is recommended herein. 

Also in Section 3.4, an equation (3-11) for evaluating proposed by 

Barr, et aT. [2-10] is given. Evaluation of fi^om this equation, howbver, 
is very dubious for thunderstorms and is not recommended. For example, a 
very large and uncommonly observed value of u* is 2 m s“^ or greater (see 
Section 3.4) which from Equation 3-11 gives a value of 0 ^^ = 2.6 m s"^. This 
value which would occur only infrequently in a normal boundary layer is a 
factor of 1.8 lower than the lowest value of measured by Steiner and 
Rhyne [2-6] bearing in mind, of course, that their data were measured at 
altitudes between 12 and 8 km. 


It is assumed, however, that Steiner and Rhyne's data will extend to 
lower altitudes and the procedure recommended for evaluating is to select 
a value in the range of 10.20 > <^wz ^ predict and owy. 

Equation 2-5 should be used. 

The length scale of the turbulence may be selected from Figure 2-5 taken 
from Reference [2-10] where 
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1.2 



Figure 2-5 Integral Length Scale Description [2-10] 
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L = L = ^ ^ 

'"'"'v ~ pOO m, z 300 m 

A y ^ 


, _ /z[0.177 + 0.00274 z]"^*^; z < 300 m 

■-w^ ■ \300 m; z > 300 m (2-5) 

Table 2-1 provides values of cyclic frequencies for different veloc- 
ities of and for spatial frequency of 0.4 and 0.004 radians/m, 
respectively. It is observed from this table that the tabulated values 
exceed 0.1 Hz only for the upper limit of and at velocities greater than 
10 m s"\ The importance of this observation to the simulation model of wind 
shear being developed herein is that the wind data utilized in the math 
model has been integrated over 10 sec periods and, therefore, contains 
oscillations in the low frequency range of less than 0.1 Hz. Therefore, it 
is proposed that a turbulence simulation for the current thunderstorm model 
have all frequencies below 0.1 Hz filtered out using a high pass filter (see 
Reference [2-12]). In many ways this result is advantageous because the 
turbulence spectra are better known at high frequencies and can be simulated 
very well; whereas, low frequency turbulence spectra are not yet well 
defined. The high frequency range corresponds to the well established 
inertial subrange [2-13] for which <j>(Q) obeys the -5/3 power law. A high 
pass digital filter program was used with the computer program given in 
Appendix 2C but is not included in the write-up because it is a staridard 
program given in Reference [2-14] and is too long to be included in this 
report. 

2.5 Statistical Model of Turbulence 

In an attempt to make some estimate of the most extreme wind shear 
which could be encountered during the expected life cycle of a thunderstorm, 
the 20 thunderstorm cases were taken as a sample for statistical purposes. 
From this sample, wind speeds along the 20 flight paths illustrated in 
Appendix 2B were statistically analyzed. Table 2-1 shows the ensemble 
average of the wind speeds along the flight paths plus the one point 
standard deviation of each wind speed and the correlations betv/een the wind 
speed components. These values are defined as shov/n below for the 
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TABLE 2-1 

ENSEMBLE AVERAGES, STANDARD DEVIATIONS AND CORRELATIONS 

FOR W^, AND 

X jr 2* 


"o 


<W > 

y 

<W^> 

0 

9.0700 

-0.6850 

0.0000 

50 

9.189b 

-1.0350 

-0.0519 

100 

9,7518 

-1.2028 

-0.1174 

150 

11.4401 

-1.3898 

-0.2904 

200 

12.3359 

-1.2510 

-0.2539 

250 

12.5787 

-0.1779 

-0.0397 

300 

11.8227 

0.2164 

-0.0670 

350 

8.8224 

1.1370 

0.6620 

400 

4.3032 

2.1131 

1.2981 

450 

1 .**57 0 

2.7131 

0.8968 

500 

-0.1681 

2.*9205 

0.4902 

^0 

"“x 


“«z 

0 

4.0341 

3.7439 

0.0000 

50 

3.4693 

3.4720 

0.1404 

100 

3.7829 

. 3.6212 

0.2641 

150 

3.9754 

3.7328 

0.4325 

200 

4.0710 

3.7115 

0,6148 

250 

3.8818 

4,9743 

1.0989 

300 

3.8244 

5,1184 

0.8581 

350 

5.4606 

4.7680 

0.9819 

400 

5,7848 

5.8223 

1.2641 

450 

5.5734 

7.3732 

0.6563 

500 

5.3830 

8.2354 

0.4774 

^0 

Vy 

Vz 

p 

Vz 

0 

3.4760 

0.0000 

0.0000 

50 

1 .9669 

0.0184 

-0.0567 

100 

1 .5427 

-0.1086 

0.1101 

150 

0.3304 

-0.7416 

-0.0176 

200 

0.1972 

-0.8220 

1 .2121 

250 

0,9971 

-0.8448 

3.8696 

300 

0.0995 

-1.4322 

2.1130 

350 

2.6395 

-2.6172 

1.2844 

400 

-4.8549 

-2.5340 

1.9091 

450 

-4.6407 

-0.5253 

-0,5719 

500 

-2.0084 

-0.1117 

-1.8076 
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longitudinal wind speed. 


■ 1 1 \ 


{2-7) 




N 

I(W, 

i=l ^i 


<Wx>r 


N - 1 


1/2 


( 2 - 8 ) 


and 




1 

^ i=l M 


<Wx>)(Wy^ - <Wy>) 


(2-9) 


The relationship for the lateral and vertical wind speed and other two 
correlations are similarly defined. 

The results of Table 2-1 allow a tri variate probability density function 
for the wind components at each point along the flight path to be con- 
structed, A trivariate distribution as desqribed in Reference [2-15] is 
given by 


p(w^, Wy- “z> ° (2-10) 

where 

C = 2TT"^/f|s..|^/2 (2-11) 

' J 


and 

- 2S,,.Z| + 2Zs,jZ,Zj; 1 = 1,2. 3 


( 2 - 12 ) 


The vertical bars designate the determinant of the matrix s^-j where 
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( 2 - 13 ) 



s 


IJ 




• • 

'^W 


and z. is the standardized variate 
= (W^. - <WT>)/s.j; j = 1.2,3 


From this function the probability of the wind vectors at each point can be 
estimated, however, the mathematics involved is too complicated to present in 
detail in this report and the reader Is referred to Reference [2-15]. 

There is considerable encouragement that a probability model of 
thunderstorms could be constructed since the components of [s^j] and the 
ensemble averages for the 20 storms are reasonably consistent functions of 
the spatial coordinate x,z. This is demonstrated in Figure 2-6 v/hich shows 
the variation of <W >, and Pw v/ along horizontal lines through the storm 
at equal increments of At = 1 min or Ax/W . The lines are at the 400, 300 

A 

and 200 m level. The data show well defined variation with time and con- 
sistent trends with altitude. Equally consistent results are obtained with 
the remaining statistical parameters. 

These well behaved statistical parameters suggest that a study to 
develop a complete statistical model of a thunderstorm be carried out. The 
results of this study would allow simulation of thunderstorms, similar to 
turbulence simulations, to be performed. Thus, a simulator could be pro- 
grammed to simulate approaches and landings through a random selection of 
thunderstorms. Moreover, the statistical model would allow an estimate of 
the most extreme wind shear associated with a thunderstorm and its frequency 
of occurrence. 

Until such a model is available, the best estimate of the probability 
of an extreme value in the thunderstorm can be made from the tabulated 
values of standard deviation (Table 2-1) and the assumption of a Gaussian 
distribution. For example, one can estimate that at the point z = 250 m on 
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the; 3® glide slope will on the average be 12.6 m s"\ there is however a 

^ -1 

32% chance that it will be ±3.88 m s of that value; a 5% chance that it 
will be ±7.76 m s“^ of that value; etc. This approach assumes that the wind 
fluctuations at neighboring points are statistically independent which is 
obviously not the case in view of the significant value of the correlation 
coefficients also shown in Table 2-1. 

Considerably more work is required to produce a satisfactory statistical 
model of a thunderstorm wind field, however, preliminary analyses conducted 
on the sample of 20 thunderstorms suggest a meaningful model is feasible and 
within reach. 

2.6 Application of the Thunderstorm Wind Shear Model 

The user of the thunderstorm wind shear model simply selects which 
thunderstorm case is of interest to his simulation problem and generates 
wind speed profiles for given values of x and z throughout the wind field. 
The selection of the thunderstorm case must consider two factors; one is the 
length of horizontal fetch over which the simulation is to be carried out 
and the other is the severity of the storm which it is desired to simulate. 

All profiles have the same vertical extent which is 500 m. However, 
they have different lengths depending upon the mean motion of the thunder- 
storm. If a simulation is to be carried out over a significantly long 
distance in the horizontal direction, one must select the storm which will 
span the distance of the proposed simulation. Table 2-2 lists the length of 
each thunderstorm based on its case number. This table can be used to 
select which thunderstorm case to use in carrying out the simulation where 
the length of horizontal extent is critical. 

The severity of the storms also differ for a number of reasons. One 
is that the maturity of the storms as they passed over the tower were in 
different stages of development, another is that the center of the thunder- 
storm may not have passed over the tower and only the fringes of the storm 
were recorded. In selecting a storm to have a severity appropriate for the 
simulation to be carried out, the user may either inspect the graphical 
illustrations shown in Appendix 2B and thus pick what would appear to be the 


24 



TABLE 2-2 


HORIZONTAL LENGTH OF EACH STORM RECORD 


Case 

Length (km) 

, Case 

Length (km) 

1 

4.00 

11 

8.16 

2 

3.28 

12 

3,64 

3 

5,71 

13 

6.35 

4 

7.69 

14 

5,13 

5 

11.43 

15 

7.41 

6 

7.27 

16 

4.44 

7 

7.84 

17 

11.43 

8 

5.56 

18 

16.00 

9 

7.41 

19 

7.14 

10 

8.51 

20 

13,33 


worst storm relative to the simulation he wishes to perform. Note, however, 
that the wind profiles shown are based on a specified flight path which 
terminates at the bottom left-hand corner of each data set and which has a 
constant 3° glide slope. Therefore, the flight path for which the winds are 
illustrated may not pass through the worst part of the storm. In this case, 
the user may examine the tables provided in Appendix 2A and by examining the 
overall wind fields select from the tables the thunderstorm case which gives 
either the largest longitudinal wind, largest vertical wind, or the largest 
lateral wind, which ever may be of interest. 

Having chosen the thunderstorm of interest, one would normally wish to 
predict what the probability of encountering such a thunderstorm would be 
and with what frequency would such storms occur. Statistical models which 
would allow this type of risk of exceedance estimate to be made are not 
available and require further research for their development. The best 
procedure to achieve some estimate of the probability of encountering such a 
storm is to utilize the information provided in Section 2.5. Table 2-1 in 
this section gives the ensemble averaged wind speeds for all 20 thunderstorm 
cases and the standard deviations of these wind speeds about this mean. 
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Again, these results are compiled for only those flight paths as specified 
earlier. The user can, by comparing his selected storm with the ensemble 
averaged value, determine the number of standard deviations his storm departs 
from the mean and in this way estimate the probability of the wind field 
magnitude. The addition or subtraction of standard deviations about the 
mean at each point along the wind field provides only a crude estimate of 
the statistics of the thunderstorm wind fields, however. This point is well 
illustrated by inspection of the coherence between wind speed components 
provided in Table 2-1. There is a very strong coherence between wind speed 
components and therefore each point in the wind field does not behave inde- 
pendently. Therefore, it is incorrect to simply add standard deviation at 
each point but in lieu of a better approach this method can be followed. 

The foregoing arguments clearly indicate that there is a need to carry 
out a more detailed analysis of thunderstorms such that a statistical model 
which would provide the extreme magnitudes of wind speed and of the wind 
shear expected to occur in a thunderstorm will be available. The frequency 
of occurrence of the extreme is also greatly needed. 

The turbulence model developed in conjunction with the thunderstorm 
wind shear model should definitely be used in carrying out any simulation 
process. The reason for this is that the data utilized to develop the wind 
shear model predict at most downdrafts of 3 m s"^ . Values of downdrafts as 
high as 15.5 m s~\ however, are reported in [2-16 and 2-17]. These values, 
however, are undoubtedly averaged over much shorter periods of time than 10 
sec for which the data presented herein are averaged. Neither reference 
gives any information on the averaging time utilized in arriving at the 
quoted value of 15.5 m s~^. Two models of thunderstorm wind fields that 
have been developed by Keenan [2-17] are tabulated in Tables 2-3 and 2-4. 

These wind fields were reconstructed from the flight data recorder of aircraft 
involved in accidents resulting from flight through severe thunderstorms. 
Inspection of these tabulated results illustrate that much more extreme 
downdrafts or downbursts as defined in Reference [2-16] occur in these data 
sets than in those tabulated in Appendix 2A. 

There is at this time conflicting data and opinions as to the maximum 
magnitude of the downdraft that can occur in a thunderstorm. Although the 
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Table 2-3 Thunderstorm Wind Field, BIO, Similar to Philadelphia/Allegheny 
Profile [2-1^; HX - Horizontal Station (ft,). Velocities Are 
Given in Knots and Height in Feet, 
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Table 2-4 Thunderstorm Wind Field, Bll, Similiar to Kennedy/ Eastern 
Profile, [2-l2; HX - Horizontal Station (ft,). Velocities 
Are Given in Knots and Height in Feet, 
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10 sec average data of Goff [2-1] will have lower values than the peak 
downdraft wind speeds reported by Fujita [2-16], the discrepancy in the 
values cannot be completely contributed to averaging time. A recent report 
by Alexander [2-18] gives a statistical summary of vertical wind speed data 
recorded at NASA's 150 m ground wind tower facility, Kennedy Space Center, 
Florida. One year of continuous around-the-clock vertical wind speed measure- 
ments were processed to determine the intensity, frequency, time of occur- 
rence, etc., of the daily maximum vertical gust. Both updrafts and downdrafts 
were studied. These values represent 0.1 sec averages and the maximum 
vertical downdraft recorded is 9.3 m s”^ although data recorded specifically 
during the hurricane Agnes indicated a downdraft in excess of 11.9 m s"^. 

Sinclair [2-19] indicates that downdrafts at 1000 m for an Oklahoma 
thunderstorm may be considerably in excess of the 15.5 m s"^ recorded in 
Reference [2-16]. Sinclair has experienced and measured downdrafts as high 
as 28 m s"^ based on a 1/25 sec averaging period. Finally, the numerical 
model of Williams, et al. [2-20] does not predict wind speed downdrafts 
greater than 10 m s"\ 

Thus, it is evident that research is needed to brin^ together the data 
currently available and to resolve the magnitude of the maximum downdraft 
which can occur within a thunderstorm. This would allow the current simula- 
tion model to be updated by superimposing turbulence fluctuations of 
realistic magnitude on quasi-steady wind fields. For the time being, however, 
to provide an estimate of wind fields which would be consistent with the 
higher values of vertical wind speed reported in References [2-16 and 2-17] 
consider the following. 

If the wind shear model based on the thunderstorm data from Goff [2-1] 

incorporates into the turbulence simulation fluctuations which are based on 

4.75 ^ s 10.2 m s"^ recorded at 12 to 8 km as described earlier in Section 

2.4, very high downdrafts will be computed. For example, taking the average 

value of = 7.5 m s"^ and adding to that the 10 sec average wind speed of 
^ -1 -1 

approximately 2.5 ms , a value of 10 m s is obtained for one standard 
deviation and a value of 17.5 m s"^ for two standard deviations. The 
reported value of 15.5 m s"^ mentioned earlier is slightly less than two 
standard deviations about the 10 sec mean. Thus, statistically it is clear 
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that downdrafts of 15 m s"^ or greater can readily occur in thunderstorms if 
the turbulence intensities at altitudes of 12 to 8 km extend to the ground. 
However, turbulence intensity is normally attenuated near the ground and it 
is not confirmed that such high values exist there. Until this is experi- 
mentally resolved, it is recommended that the model of thunderstorm wind 
shear based on the extensive data from Goff [2-1] should have turbulent 
fluctuations superimposed with a standard deviation of 7.5 m s~V. The 

downdraft magnitudes reported in accident investigation will then be achieved 
in the simulation. 

It is anticipated that the proposed turbulence simulation model based on 
the work of Reeves, et al. [2-3] will provide a realistic turbulence simu- 
lation to accompany the wind shear model proposed in this report. However, a 
better simulation of the large downdraft fluctuations that occur in a thunder 
storm can be achieved by utilizing a model of turbulence which includes 
coherence between different levels of the atmosphere or between different 
positions in the storm. Such a model has been preliminarily developed by 
Perlmutter and Frost [2-21]; however, this model requires further perfection 
and the coherence function associated with thunderstorms must be developed to 
permit its use. 

An alternative explanation of why the data of Goff does not contain 
downbursts of the intensity reported in [2-16 and 2-17] is that the 20 
thunderstorms investigated may not contain a "spearhead echo" type storm 
[2-16] or the measured data may not encompass the downdraft portion of the 
storm. Fujita [2-16] defines a downburst as = 3.6 m s“^ at the 300 m 
level. The averaging time related to this wind speed is not specified, and 
it is not known whether this represents a peak gust or a value averaged over 
some interval of time which is undoubtedly less than 10 sec. The data of 
Goff [2-1] do not indicate any values of that equal or exceed the 3.6 m s' 

definition of a downburst. There are a few values, however, that do approach 
the 3.6 m s'^ level; for example, thunderstorm cases 8 and 11 shown in Table 
2-5 which lists the maximum and minimum values of recorded at the 300 m 

level in the data utilized to construct the thunderstorm wind shear model 
given in this report. The fact that no downbursts are recorded is not sur- 
prising, however, because obviously the chance of a downburst being directly 

over the tower the instant of maximum intensity is extremely small. 

34 



TABLE 2-5 

MAXIMUM AND MINIMUM VERTICAL VELOCITIES AT 300 tn 
LEVEL IN DATA OF GOFF [2-1] 


Case 

W (m s~^) 

^max 

W (m s”^) 

^min 

Case 

W_ (m s'^) 
max 

W_ (m s"‘) 
^mi n 

1 

1.6 

-2.5 

11 

1.8 

-3.0 

2 

3.0 

-1.2 

12 

2.0 

-0.9 

3 

2.6 

-1.1 

13 

3.2 

-0.7 

4 

3.0 

0.7 

14 

2.1 

-1.0 

5 

3.0 

-1.3 

15 

2.0 

-0.5 

6 

3.7 

-1 .5 

16 

2.0 

-1.0 

7 

1.2 

-2.1 

17 

1.1 

-0.5 

8 

4.0 

-3.0 

18 

1 .0 

-1.0 

9 

4.1 

-0.1 

19 

2.0 

-2.0 

10 

3.0 

-1.2 

20 

0.2 

-1.1 


Fujita [2-22] points out the probability of an airport being under the 
influence of a spearhead echo is very low and is probably less than 2 percent 
of the thunderstorm probability. Moreover, he notes that the location of 
aviation hazards for the extreme downburst is limited to only a fraction of 
the spearhead echo area. With this in mind it is relatively easy to believe 
that the data of Goff do contain a few storms of the magnitude approaching 
those types defined as spearhead echo by Fujita. Moreover, since several of 
the storms clearly record downdrafts approaching the magnitude of downbursts 
for 10 sec averaged wind speed, it is envisioned that significantly higher 
downward gust for a shorter, say, 3 sec average are contained in the original 
data. 

The data set illustrated in Tables 2-3 and 2-4 is being punched on cards 
and will be included with thunderstorm data now stored on magnetic tape. It 
is believed, however, that the current data set utilizing the appropriate 
turbulent intensities will provide a valid simulation of thunderstorms for 
flight/hazard simulation studies. The larger selection of thunderstorms, 20 
in number, enables a manned flight simulation study to provide the pilot with 
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several different thunderstorm situations that do not duplicate those that he 
has negotiated in previous tests. Thus, he cannot learn a given storm. By 
selecting the same sequence of thunderstorms, a second pilot involved in the 
same test program can be exposed to the identical test pattern utilized for 
the first pilot with neither pilot flying the same storm more than once. 

Finally, it is called to the attention of the user, that in carrying out 
a simulation where avionics using ground wind speeds as inputs are being 
studied, the storm must be considered to be passing over the airport or 
anemometer from which the ground wind is being determined at the speed of the 
gust front, W . Therefore* the length of the thunderstorm record used in the 

X 

simulation must be sufficiently long that the location of the assumed 
anemometer at the airport has not moved out of the range of data set during 
the time taken by the aircraft to complete its approach. 
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SECTION 3.0 


NEUTRAL AND STABLE BOUNDARY LAYERS 

The data set utilized in formulating a mathematical model for wind 
shear or wind fields associated with neutral and stable boundary layers 
relies mainly on the work of Clarke and Hess [3-1]. Although numerous 
analytical models for boundary layers, both undOr neutral and stable condi- 
tions w4re available in the literature, these were in general based on the 
assumption of constant stress which is valid only to the first 50 or 100 m 
of the atmospheric boundary layer. Other analytical models which did not 
evoke this assumption and included the influence of turbulence of the atmo- 
spheric boundary layer were generally highly mathematical and required 
numerical solution with a computer. Therefore, the data set of [3-1], which 
are presented in the form of contour maps of constant longitudinal and 
lateral wind speeds as a function of height, z, and the stability of the 
atmosphere expressed by the stability parameter, u» allowed tabulation of 
the wind fields on a gfid system. These tabulated data were then coupled 
with a table lookup computer program to provide the fast time wind speed 
model for flight hazard simulation. 

None of the data sets or mathematical models available in the literature 
for boundary layers (see Reference [3-2]) incorporate horizontal variation 
of wind speed. Hence, the wind shear models for neutral and stable boundary 
layers reported herein depend only on the vertical scale, z. It should be 
borne in mind, however, that terrain features indigenous to each landing 
site can influence the spatial gradient of the wind speed horizontally and 
also introduce a possible vertical wind speed component. Fortunately, most 
airport terrain is relatively flat and not surrounded by high mountains or 
sharp cliffs which create flow disturbances, both vertically and horizontally 
in the wind. Therefore, the wind shear model described in the following 
section is expected to present a valid simulation of wind shear in terminal 
areas of the typical airport. However, indiscriminate application of the 
model is cautioned against until further studies of wind fields over irregular 
terrain become available. 
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3.1 Data Source 

Boundary layer data are presented in the form of horizontal and lateral 
wind vector components. The boundary layers are described over flat uniform 
homogeneous terrain and therefore have no vertical component of wind speed 
or dependence on the x-coordinate. Of course, when turbulence fluctuations 
are superimposed as described in Section 3.4, a vertical component can 
occur. The wind speed profile is thus dependent only on the height, z, and 
on the stability of the atmosphere characterized by the stability parameter, 
p. The data set used in the boundary layer model is based on the results of 
the extensive wind speed measurements reported in Reference [3-1]. Also, 
the influence of baroclinicity is neglected and is justified for this data 
in Reference [3-1]. 

Wind speeds and temperatures were measured over a very flat, smooth 
surface having a mean surface roughness parameter, z^, of approximately 3.5 
mm. Hourly double theodolite observations of pilot balloons were taken at 
five stations for 40 days, Micrometeorological observations were taken at 
two of the five stations. Wind profiles were measured at 0.5, 1, 2, 4, 8 
and 16 m while temperature differences were measured at one station between 
1 and 2 m, 2 and 4 m. Values of u* were estimated by the drag coefficient 
method and values of surface heat flux by means of temperature and wind 
profiles. Surface pressures were measured at each of the five stations and 
radiosonde measurements were made every three hours at the central station. 

In addition, three hourly surface pressure and temperature data from 14 
stations at distances up to 350 km away were used to augment the data obtained 
by the research expedition. 

Geos trophic winds were determined from the data by using second order 
curve fitting procedures on the pressure data. Thermal winds were determined 
using temperatures from a surface network reported every three hours and 
from a temperature field measured each day at 1500 hours. Interpolation 
from the three-hour data set were primarily used in assessing baroclinicity 
which was shown to be small. A second set of thermal wind estimates based 
on twice daily radiosonde data from a network of five stations were also 
made. Reference [3-1] reports that surface geostrophic winds were well 
determined as evidenced by the high correlation (93 percent) with observed 
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winds; however, no claim of high accuracy in estimating thermal winds is 
made in the report. 

These data thus reduced were reported by Clark and Hess [3-1] in the 

/s 

form of contour maps of winds as a function of dimensionless height, z = 
zf/u*Q, and of the stability parameter, y = icu*/fL'. The contours of the 
map are lines of constant aQ = W (z = 0.15) - W (z) and AW * W„(z) - 
Wy(z = 0.15), respectively, where W is the dimensionless wind speed W/u*. 

The symbol f is the Coriolis parameter (here, treated positive in both 
hemispheres with a right-hand coordinate system implied in the Northern 
hemisphere). For use in this report, the data were converted to a 34 x 11 
grid as illustrated in Figure 3-1. The data were then stored on magnetic 
tape and a computer look-up routine developed for interpolating the values 
of and W^ for given values of z and y. The tabulated data are given in 
Appendix 3A, illustrative wind speed profiles are given in Appendix 3B, and 
a computer program for looking up and interpolating the data is given in 
Appendix 3C. 

3.2 Presentation of Data 

/V 

The data are presented in a right-hand coordinate system with W 
positive in the direction from left to right and W positive into the plane 
of the paper. 

The values of y range from -333.34 to 216.67 in increments of Ay = 

16.67. This unusual increment size was chosen for convenience in extracting 
the data from the contour plots. Values of z range from 0.001 to 0.15 in 
increments of Az = 0.0149. 

Because of the similarity and scaling laws used in reduction of the 
data, values of wind below z = 0.001 are not given. To establish the profile 
from z = 0.001 to zero the log-linear wind speed profile was used (for a 
description of the log-linear wind speed profile see Reference [3-2]). 

Values of Q (z) below z = 0.001 are determined from 

X 

W (z) = i {£n(Ro z + 1) + 4.5 zy/K}; y > 0, 0 s z < 0.001 (3-1) 

X K 

where Ro is the Rossby number defined as Ro = u*/fz^. The variable z^ is the 
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Figure 3-1 Grid System Superimposed on the Data of Clarke and Hess 
[3-1] 
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empiricany determined surface roughness. Typical values of are given in 
Reference [3-2]. The value of is zero for 0 < z < 0.001. Introducing 

z - O.bOl into the equation gives the value of wind speed to which all the 
tabulated longitudinal wind speeds are referenced. 

3.3 Mean Wind Speeds 

This section considers values of the mean wind fields computed from the 
mathematical model of the neutral and the stable boundary layers. Although 
the data include the range of unstable condition, i.e., -334.34 i y <0, 
strong wind shear is not normally associated with unstable boundary layers 
and hence no values of winds for this range of y are given. 

The wind data are averaged over a 10 minute period or longer and thus 
represent mean values. A model of turbulence for the neutral and stable 
boundary layer is given in Section 3.4 which is superimposed to give the 
random instantaneous wind speed. 

The grid system used in storing the data is numbered from the bottom 
left-hand corner. The numbers increase from left to right in the direction 
of increasing y and from bottom to top in the direction of increasing z. 

Thus column 1 corresponds to y = -333.34 and column 34 corresponds to j 
y = 216.67 where row 1 corresponds to z = 0.001 and row 11 corresponds^ to 
z = 0.15. 

3.3.1 Tables of Wind Speed 

Tables of the longitudinal and lateral wind speeds W (z) and W (z) for 

X y 

neutral and stable boundary layers, i.e., y ^ 0, are given in Appendix 3A. 
These values are computed from the tabulated wind differences from the 
relationship 

Qx(z) = - AWj^(z) + Sjz = 0.001) (3-2) 

and 

Wy(z) = AWy(z) - AQy(z = 0.001) (3-3) 
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where from Equation 3-1 evaluated at z = O.OOT, 

W (z = 0.001) = x"^[Jin(0.001 Ro + 1) + 0.01125 y] (3-4) 

A 

3.3.2 Illustrations 

The longitudinal and lateral wind speeds which would be encountered 
along a 3® glide slope are plotted for various stability conditions and 
Rossby numbers in Appendix 3B. Height is expressed in meters, m, on the 
vertical scale and wind speeds are expressed in meters per second, m s , on 
the horizontal axis. 


3.3.3 Computer Program 


A computer program has been written which computes with inputs of height, 
z, and stability parameter, y, the longitudinal, W^, lateral, W^, wind speeds 
and the wind gradients 9W^/3z and 3Wy/3z. The computer program also requires 
as input the friction velocity, u*, the Coriolis parameter, f, and surface 
roughness, z^. All velocities are input and output in m s"^ and lengths in 
meters, f, is introduced in s"^. 


This computer code can be used as a direct subroutine input to existing 
computer programs for flight simulators or computer programs of airplane 
motion in variable wind fields. The option of superimposing turbulence on 
the mean wind speed is available by appropriate specification of control 
variables. A complete description of the computer program is given in 
Appendix 3C. 


3.4 Turbulence Model 

A turbulence model has been developed for use with the neutral and 
stable boundary layer data which employs turbulence kinetic energy spectra 
developed by Kaimal [3-3], as shown in Figure 3-2. Kaimal gives the following 
functional form for the spectra. 


n(j)^{n) 0.164 n/riQ 

al 1 + 0.164(ti/i1q)^^^ 


(3-5) 
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4> (n) Power spectral density distribution in a direction, s 

^ - 1 
n Cyclic frequency, s 

0 ^ Turbulence intensity or root-mean-square fluctuation in 

a direction, m s 

n Reduced frequency nz/W (z) 

A 

z Height above ground, m 

Mean wind speed at the height z, m s"^ 

a Designates component of velocity fluctuation (either W , w 

X y 

or w^) 

The influence of atmospheric stability enters through the variable 
which is a characteristic reduced frequency and is a function of Richardson's 
number as shown in Figure 3-3. These values of are for the range 
0.04 < Ri < 0.20. 

Values of for the neutral boundary layer where Ri = 0 are recommended 
in Reference [3-4] as 

(n„)„ = 0.01444; = 0.0265; (n„)„ = 0.0962 (3-6) 

X y z 


Thus for the complete range of 0 < Ri < 0.2 the approximate relationships are 
used in this report are 



'0.0144 

0 

o 

X 

II 

0.5 Ri 

0.029 

1 

'0.0265 

0 

It 

o 

1.5 Ri 

0.018 


'0.0962 

0 


.2.8 Ri 

0.035 

The relati 

i onshi p 

between Ri 

Ri = (zy/i 

<)[1 + ' 

4.5 zp/k]"^ 


< Ri < 0.029 
s Ri < 0,2 

< Ri < 0.018 
^ Ri < 0.2 

< Ri < 0.035 

< Ri < 0.2 

and y is 


(3-7) 

(3-8) 

(3-9) 


(3-10) 
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Figure 3-3 Correlation of tiq with Richardson Number for Stable 
Boundary Layers 





The turbulence spectra data of Kalmal [3-3] given above appear to be the most 
comprehensive presentation of atmospheric spectra associated with the stable 
boundary layer and are taken herein to represent the state of the art. 

To utilize Equation 3-5 a value of turbulence intensity must be deter- 
mined. Figure 3-4, taken from Reference [3-5], gives the turbulence intensity 
of the vertical wind speed component, , nondimensional ized with u*, plotted 

A 

as a function of nondimensional altitude, zy/ic. For neutral conditions, 
zy/K = 0 and the ratio of vertical turbulence intensity, to the friction 
velocity, u*, becomes 1.3. 

At zy/K = 1.22, the turbulence intensity vanishes as the atmospheric 
boundary layer becomes so stable that essentially laminar flow is achieved. 

In general, decreases with increasing stability. 

A curve fit of Figure 3-4 between 0 < zy/K < 1.22 is 

1.3 - 0.13(zy/K)°-^ 0^ zy/K <1.0 

6.49 - 5.32(zy/K) 1 .0 < zy/K < 1 .22 



The value of can thus be determined for given values of u*, y and z. 

No satisfactory mathematical description of how 0^ and vary with 

large scale surface features nor how they vary with atmospheric stability is 
available. Barr, et al. [3-5] propose that the ratio /Oy^ be treated as a 
function of altitude according to the following 



/a 


w 


z 


[0.177 + 0.832 z/z.] 

1.0 


z < z. 


z > z^. 


(3-12) 


z^ = 300 m or 1000 ft. Reference [3-5] also proposes that 

This relationship does not result in satisfactory agreement with Equation 

3-5, however, and the relationship 



/a 


w 


z 


[0.583 + 0.417 z/z^.]"®*® 

1,0 


z < z. 


z > z. 


(3-13) 
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Figure 3-4 Vertical Turbulence Intensity 



, is proposed by Frost [3-4], Equation 3-13 is developed identical to Equation 

3-12 but assumes that is less than near the ground. The assumption 

y X 

is that for neutral conditions ” 2. 6/2. 0/1. 3 which is consistent 

with a number of reported results. 

Values of and are plotted in Figure 3-5 as a function of 

dimensionless height to facilitate computation of these values. Inspection 
of Equations 3-11, 3-12 and 3-13 show that for the neutral atmosphere 

a = 2.6 u*, a = 2.0 u* and a = 1.3 u* (3-14) 

X y z 

Figure 3-6 shows longitudinal turbulence spectra for varying degrees of 
stability. Turbulence in the atmosphere is observed to decrease with 
increasing stability. This is most evident at the low frequency range and is 
to be expected since the thermal effects which create stable boundary layers 
depress large scale turbulent motion. These spectra, however, are reported 
to be influenced by surface terrain features, particularly in the longi- 
tudinal and lateral components. No mathematical models which account for 
terrain effects are available. The spectra having been measured for rela- 
tively homogeneous terrain which is characteristic of the majority of 
airports should give a valid representation of most terminal areas. For air- 
ports located near unusual terrain features such as mountains or cliffs, 
fluctuations in the longitudinal and lateral components of the wind may, 
however, be higher than simulated by the proposed turbulence model . 

The spectra represented by Equation 3-5 do not have a rational form and, 
consequently, are difficult to use in turbulence simulation schemes [3-6]. 

To overcome this difficulty, the spectral data in Figure 3-2 were adjusted to 
fit a modified Dryden spectrum having the functional form 

ni(n) _ (3_15) 

0^ 1 + C^Cn/Ho)^ 

The constants C-j and are adjusted such that the modified Dryden curve fits 
Kaimal 's data [3-3]. Values of C-j = 0.1580 and C2 = 0.0694 were determined 
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Figure 3-5 Longitudinal and Lateral Turbulence Intensity 
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Figure 3-6 Influence 


to give the best curve fit. Thus, the power spectral density function chosen 
for the turbulence simulation model 


nt(n) _ 

1 + 0.0694(n/r,j,)^ 


(3-16) 


Figure 3-7 compares the modified Dryden spectra with Equation 3-16. Although 
this curve fit does not provide a good representation of the higher frequency 
spectra, it does fit the data quite well in the lower frequency range which 
is expected to have the most significant influence on the flight of aircraft. 

This turbulence model utilizing the turbulence spectra given by Equation 
3-16 and the z transformation technique (see Neuman and Foster [3-7]) has 
been developed and can be coupled with the mean wind fields for the stable 
and neutral boundary layers to give a random fluctuating field. The model, 
however, is linear and results in a Gaussian distribution of the wind speeds 
in the atmosphere which introduces a small error into the simulation (see 
Reeves, et al . [3-6]). 

To illustrate the influence of turbulent fluctuations on the velocity 
profile experienced by an aircraft during landing in a stable boundary layer. 
Figures 3-8 through 3-13 have been prepared. These figures illustrate stable 
boundary layer wind speed profiles seen by an aircraft on a 3° glide slope 
with turbulence, computed by the simulation technique, superimposed. The 
inertial aircraft velocity is 64 m s“^ which corresponds to a sink rate of 
3.35 m s“^. The time increment used in the turbulence simulation was taken 
as 0.15 sec which results in a turbulent fluctuation being superimposed at Az 
increments of 0.35 m along the flight path. The figures show some interest- 
ing results relative to the stability of the boundary layer. In these 
figures, u* has been held constant and p is increased in sequential order. 

The first velocity profile (Figure 3-8) corresponds to the neutrally stable 
layer where y is essentially zero. One observes that the aircraft encounters 
turbulence from the 500 m level to the ground and that the turning of the 
boundary layer is reasonably small after 100 m. The second figure (Figure 
3-9) illustrates a somewhat higher level of stability, y = 25. The longi- 
tudinal wind speed is larger because, as noted earlier, the computed 
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3-7 Comparison of the Modified Dryden Spectrum with Kaimal's Emperical Gurve Fit [3-3] 





450o0 



Figure 3-8 Near NeutreL'l Boundary Layer with Turbulence Superimposed 





Figure 3-9 Stable Boundary Layer with Turbulence Superimposed 



- 500,0 



Figure 3-10 Stable Boundary Layer with Turbulence Superimposed 






Figure 3-12 Stable Boundary Layer with Turbulence Superimposed 



- 450,0 



Figure 3-13 Stable Boundary Layer with Turbulence Superimposed 


profiles are based on the same friction velocity, u*, and surface roughness, 
Zq. The interesting feature of this figure is that at large values of z, 
i.e., where z/L' is large, the mechanical turbulence is damped out by buoyancy 
induced turbulence, however, as z/L' becomes small due to the aircraft 
approaching the ground, the situation is reached where the atmospheric 
boundary layer returns to a neutral condition and mechanical tiirbulence 
dominates the buoyancy induced turbulence. In Figure 3-9 this occurs at 
approximately 125 m., The airplane thus flies from a region of rather quiet 
flow to a sudden and rapidly increasing turbulent intensity. The turbulence 
intensity becomes larger as y increases, however, this is not due to the 
effects of stability but due to higher mean velocities resulting from a 
constant u*. For consecutive figures, as the stability of the boundary layer 
increases one observes that the transition from the essentially laminar flow 
to the highly turbulent region approaches lower and lower levels. It should 
also be noted from this sequence of figures that the strongest directional 
shear occurs at intermediate values of y, whereas, the strongest linear shear 
occurs at the larger values of y. 

An alternate way of computing the profiles would be to assign the same 
wind speed at a certain level and adjust u* accordingly. Physically this 
corresponds to conditions encountered at an airport where the wind speed 
measured with an anemometer at a fixed height, say, at the TO m level would 
record identical values but the wind shear would be appreciably different due 
to the associated stability conditions. Figures 3-14 through 3-19 show com- 
puted profiles having a common value of 10 m s"^ wind speed at the 10 m 
level. The value of u* used in the computation is related to the stability 
parameter, y, by 

u* = 0.58 - 0.0016 y 

The same characteristics of the wind speed profiles are observed, how- 
ever, the higher wind speeds now occur at intermediate values of y rather 
than at the higher values as in the sequence of Figures 3-8 through 3-13. In 
turn, one observes that the directional shear is also largest at intermediate 
values of y which is in correspondence with the former sequence of figures. 
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Figure 3-14 Near Neutral Boundary Layer with Turbulence Superimposed 



■^oOoO-i / u*=a. 56 m 



Figure 3-15.. Stable Boundary Layer with Turbulence Superimposed 




Figure 3-16 Stable Boundary Layer with Turbulence Superimposed 
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Figure 3-18 Stable Boundary Layer with TurbulenGe Superimposed 
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Figure 3-19 Stable Boundary Layer with Turbulence Superiniposed 



3.5 Statistical Model Relative to Risk of Exceedance 


The wind shear associated with the stable and neutral boundary layers 
are dependent upon the variables, p and u*, or the mean wind speed, W , since 
they are directly related. A statistical description of these two parameters 
will enable the user of the wind shear model to provide an estimate of the 
probability and frequency of exceeding a prescribed value of wind shear. 

In order to establish this statistical description some estimate of the 
probability of a given stability condition occurring simultaneously with a 
given value of u* is needed. An approximate analysis to achieve this goal is 
proposed as follows. From statistical theory the probability of u and u* 
occurring simultaneously is equal to the probability of y given u*, 

P(p > y /u*), times the probability of u*, P(u* > u* ), [3-8]. This is 

H r 

referred to as a conditional probability. 


P(y > bp*, u* > u*p) = P(y > yp/u*) P(u* > u*p) (3-17) 

The probability of y given a wind speed, W , which can be directly 

A 

related to u*, has been estimated by Barr, et al. [3-5], They give the 
probability of Richardson number, Ri, for a given wind speed, W^, where both 
W and Ri are evaluated at 6.1 m, shown in Figure 3-20. Note that these 

X 

curves are highly interpolated and based on only two sets of data. The 
Richardson number, Ri, is related to y by the relationship 


Ri 


y/u* 

K/zf + 4.5 y/u* 
0.18 

KU* 


0 < Ri < 0.18 
Ri > 0.18 


and W is related to u* and y by 

A 


(3-18), 



i^n 

< 


K 


In 


KU^ 


+ 4.5 






^ + 1 + 4.5 + 5.5 Jin 


+ 1 




0 ^ Ri < 0,18 

(3-19) 

Ri > 0.18 
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Figure 3-20(c) Cumulative Percent Frequency of Occurrence of Ri at Given Wind Speeds, 
12 to 15 m s”^(28 to 32 mph) 



Assigning values of z =6.1 m, = 0.01 m, f = 10"^ s"^ , k = 0.4, Equations 
3-18 and 3-19 have been plotted in Figures 3-21 and 3-22 for convenient 
reference. 

Using the above results the designer may select a value of u and a value 
of u* which gives a wind shear condition in which he is interested. The 
probability of y given u* is then obtained directly from the figure. 

The probability of u* can be determined from the distribution of mean 
wind speeds given by Frost [3-4] as interpreted from Justus, et a1 . [3-9]. 
This reference shows that the probability of a mean wind speed greater than a 
prescribed value W occurring during the year is described by a Rayleigh 
distribution, i.e.. 


P(W > Wp) 




(3-20) 


Reference [3-4] shows that this equation is a general representation of wind 
speeds measured at 138 airport sites throughout the United States where an 
average value of the scale factor, c, is 4 m s“^ ± 0.9 m s"^. Since u* is 
related to mean wind speed by Equation 3-19, the probability of u* is 
directly related to the probability of W . 

X 

With the information provided above, the user of the neutral and stable 
boundary layer wind shear model can estimate what the probability of a given 
wind shear occurring over an expected lifetime is by utilizing the prob- 
ability of y given u* from Figure 3-20 and multiplying that by the 
probability of u* given by Equation 3-19. To illustrate the use of the above 
model, consider the following example. A user wishes to determine the 
probability of y being equal to 200 given that u* is equal to 0.5 m s“^. 

From Figure 3-21 find Ri = 0.164 corresponding to y = 200 and u* = 0.5 m s~\ 
i.e., y/u* = 400. From Figure 3-22 determine W /u* = 23 and therefore W = 

— ^ X X 

11.5 m s for the given values of u* and y. From Figure 3-20 the probability 
of Ri > to 0.16 at 11.5 m s'^ is 1.0 - P(Ri > 0.16; W^ = 11.5 m s"^) = 0.12. 
From Equation 3-20 the probability of W ^ 11.5 m s“^ is P(W > 11.5 m s"^) = 
0.0003. The probability of both y > 200 and W^ > 11.5 occurring in a given 
year is from Equation 3-17 equal to 0.003%. 
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Figure 3-21 Variation of Ri with p/u* at z=6ol 




Figure 3-22 Variation of Wind Speed at 


To establish the risk of exceeding a prescribed wind shear in a given 
period of time we proceed as folTows. The probability, p, that p will be 
exceeded in any one year is related to the return period T by the relation- 
ship Tp = 1. The probability that a value less than or equal to p will occur 
in one year is q = 1 - p. In a period of N years the probability Q that a 
value less than or equal to p will occur is 

q = q^* (3-21) 

The probability or risk that a value greater than p will occur at least once 
in a period of N years is 



= 1 - (1 - p)'^ (3-22) 

Hence continuing the foregoing example, the probability that p = 200 and 
u* = 0.5 m s"^ will occur at least once in a period of 25 years is from 
Equation 3-22 equal to 3.4%. 

To give the reader some feel for the nature of the wind shear associated 
with a given probability of occurrence. Table 3-1 has been prepared. This 
table gives the risk of exceedance associated with each of the wind shear 
conditions depicted in Figures 3-8 through 3-18. 

3.6 Application of Wind Shear Models for Neutral and Stable 
Boundary Layers 

Application of the wind shear model for the neutral and stable boundary 
layer proceeds by first selecting stability conditions for which the simula- 
tion is to be carried out. The parameters required are the stability 
parameter, p, the friction velocity, u*, and the surface roughness, z^. The 
rule-of-thumb for the effect of these individual parameters is that high 
values of u*, p and all result in larger linear shear. On the other hand, 

intermediate values of p, say, on the order of 50, give the largest direc- 
tional shear within the lower 500 m of the atmospheric boundary layer. 
Inspection of Figures 3-8 and 3-9 clearly illustrate this. 
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TABLE 3-1 


RISK OF EXCEEDANCE ASSOCIATED WITH WIND SPEED PROFILES 
SHOWN IN FIGURES 3-8 THROUGH 3-19 


U* 

y 

Ri 


P(W^) 

P(Ri/Wj,) 

P(Ri & W^) 

Risk of 
Exceedance 
in 25 Yrs. 

Do 50 

200 

0.18 

13.64 

8,85x10“^ 

0.70 

6.2x1 0”^ 

0.16% 

0.50 

150 

0.17 

12.24 

8.61x10"^ 

0.74 

6.4x10"^ 

0.16% 

0.50 

100 

0.15 

10.83 

6.54x10"^ 

0.65 

4.3x10*^ 

1.07% 

0.50 

75 

0.14 

10.13 

1.64x10"^ 

0.64 

1.0x10”^ 

2.47% 

0.50 

50 

0.12 

9.43 

3.88x10"^ 

0.58 

2.3x10”^ 

5.59% 

0.50 

25 

0.080 

8.72 

8.61x10"^ 

0.56 

4.8x10"^ 

11.33% 

0.50 

1 

0.005 

8.05 

17.47x10"^ 

0.42 

7.3x10“^ 

16.74% 

0.25 

200 

0.35 

6.32 

8.96x10“^ 

0.72 

6.5x10"^ 

0.002% 

0.34 

150 

0.20 

6.21 

1.09x10"^ 

0.61 

6.6x10"^ 

0.02% 

0.42 

100 

0.16 

12.11 

1.04x10"^ 

0.68 

7.1x10"^ 

0.18% 

0.46 

75 

0.14 

11.41 

2.92x10"^ 

0.64 

1.9x10'^ 

0.47% 

0.50 

50 

0.12 

10.71 

7.72x1 0"^ 

0.62 

4.8x10“^ 

1.19% 

0.56 

12.5 

0.04 

9.65 

2.95x10"^ 

0.59 

1.7x10’^ 

4.16% 

0.58 

1 

0.0042 

9.33 

4.34x10"^ 

0.54 

2.3x10’^ 

5.59% 


Having selected the parameters of interest, one then determines the pro- 
bability that these conditions will occur and what the risk is of exceeding 
the values that are planned for the simulation. If it is desired to simulate 
severe cases of wind shear, then values of y and u* should be chosen to have 
large magnitudes which have correspondingly lower risk. Therefore, a low 
value of risk is acceptable. Simulating average daily conditions, however, 
allows a higher risk of exceeding the prescribed stability conditions to be 
accepted. The percent risk that the user of the wind shear model is willing 
to accept is subject to his own judgment based on the engineering application 
for which the simulation is being carried out and on the consequences asso- 
ciated with the test not being of sufficient severity. 

One must also bear in mind that the risk of exceeding a prescribed value 
computed in Section 3.4 and tabulated in Table 3-1 are on a per airport per 
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year basis; consequently, although the probability of u* = 0.5 and y = 200 
occurring simultaneously is 6.2 x 10’’^'*^^ (see Table 3-1), for 10,000 airports, 
this condition will occur at 0.62 of an airport or at least one airport per, 
year. In view of the fact that stability conditions may endure for one or 
two hours at an airport having heavy traffic, several pilots will be exposed 
to these conditions every year. 

Use of the turbulence simulation routine to be superimposed upon the 
steady state winds requires no additional specification of atmospheric param- 
eters, however, the time step increment of the random signal, DT, in 
subroutine STB, must be specified. In general, its value is set equal to the 
time step used in the calling program carrying out the numerical integration 
of the equations of motion for the flight dynamics of the aircraft. On the 
other hand, it can be shown [3-7] that 

T < 4 X 10"^Ti^W^/z 

but since this value varies with altitude an average value on the order of 
0.01 sec or less is generally recommended for use with the turbulence simula- 
tion routine. 

In carrying out a realistic simulation, it is recommended that the 
turbulence be used with the mean wind profile. It is evident from Figures 
3-8 through 3-19 that simulation of high stability conditions, that is, large 
y, is accompanied with a transition from a rather smooth wind aloft to 
extremely severe turbulence near the 200 to 100 m range, depending upon the 
stability. This sort of turbulence phenomenon is expected to be quite real- 
istic judging from comments with commercial airline pilots and from 
experimental data available in the literature. Super-positioning turbulence 
on the wind fields, however, does require that a computer be available for 
carrying out the computations of the wind speeds. This does not impose a 
hardship because normally, a computer is involved whether the user is per- 
forming a manned flight simulation or programming a fast time computer 
solution. 

Without turbulence, however, one can construct a wind speed profile for 
the neutral or stable boundary layer manually from the tables given in 
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Appendix 3A. As an example, determine the wind speed profile for p = 50, 
u* = 0.5 and = 0.05. To establish the profile one first computes the 
Rossby number Ro = u*/z^f = 1 x 10^. This value corresponds to Table 3A-1 in 
Appendix 3A. Selecting the column labeled y = 50 gives the wind speed 


profiles 

/s 

/s 

z 

W 


0.15 

X 

63.4 

■y 

-32.0 

0,10 

63.4 

-27.2 

0.05 

59.3 

-17.0 

0.005 

29.5 

-4.8 

0.0001 

24.1 

0.0 

Dimensional values 

/\ ^ 

are obtained from z = u*z/f, W = W u* and W^ 


hence 


z 

W 

W 


X 

y 


-1 

-1 

m 

m s 

m s 

750 

31 . 7 

-16.0 

500 

31.7 . 

-13.6 

250 

29.6 

-9.5 

25 

14.7 

-2.4 

5 

12.1 

0 


Thus, wind speed profiles can be established very simply. 
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SECTION 4,0 


FRONTAL WIND SHEAR 

A model’ of frontal wind shear has been developed in exactly the same 
manner as the thunderstorm model. Unfortunately, the number of data sets 
available for processing are only two, which are shown in their original 
format in Figures 4-1 and 4-2. These data were provided courtesy of 
Ms. Judith Stokes and Mr. Jean Lee from the National Severe Storms 
Laboratory in Norman, Oklahoma, The data represent cold frontal passage in 
the vicinity of the 500 m tower and the reduction of these data sets is the 
same as used by Goff [4-1] for thunderstorm gust front cases (see Section 
2.2). The data were measured on January 28, 1977, labeled Case 1 in this 
report, and on December 10, 1976, labeled Case 2. 

No other data which would provide the necessary detailed wind speed pro- 
files over at least a two-dimensional plane in space were found during this 
study. Additional data which is on magnetic tape at the National Severe 
Storms Laboratory have not yet been processed. No warm front wind shear data 
from which a model for flight/hazard simulation studies can be constructed 
appears to be available. 

The simulation model developed from the limited data is constructed in 
the same manner as the thunderstorm models. That is, a grid system is super- 
imposed on the contours of constant wind speed shown in Figures 4-1 and 4-2. 
In these two cases a 113 x 11 grid is used. The data were again punched on 
cards and stored on magnetic tape. The only difference between these data 
and those of the thunderstorm cases is the larger grid having 113 columns 
rather than 41 which represents a total horizontal distance of approximately 
17.2 km in Case 1 and 13.2 km in Case 2, This represents approximately 30 
minutes of recorded data. 

The following section. Section 4.1, describes the quasi-steady wind 
speed profile working data developed for the large scale frontal wind shear. 
The wind speed is referred to as quasi -steady in view of the fact that, as 
with the thunderstorm data, it has been averaged over 10 sec time intervals 
and thus contains fluctuations in the wind speed of 0.1 Hz or smaller. It is 
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Figure 4>1 Wind Speed for Synoptic Front Cold Air Out-flow 









Figure 4-2 Wind Speed for Synoptic Front Cold Air Out-flow 















proposed that the same turbulence simulation model for thunderstorms be used 
for the frontal wind shear. Since absolutely no data was found which would 
give information pertaining to turbulence intensities and characteristic 
turbulent length scales in a major frontal flow, the lower vaTues associated 
with thunderstorms are recommended as input to the model. 

4.1 Quasi -Steady Wind Speed Profile Grid System 

The wind fields shown in Figures 4-1 and 4-2 were fit to a 113 x 11 
point grid system. Wind speeds at each grid point were tabulated and punched 
onto computer cards. The data were later stored on magnetic tape. 

The grid system is numbered from the left-hand bottom corner. The 
numbers increase from left to right in the positive x-direction and from 
bottom to top in the positive z-di recti on of the original data. The wind 
speeds are given in units of meters per second, m s~^, with W being positive 
in the direction of frontal motion, W^ being positive upward and being 
positive into the plane of the paper. The frontal speeds of the tv/o storms 
are 9.3 m s*^^ and 7.1 m s"\ respectively. 

The wind speeds, W , shown in Figures 4-1 and 4-2 are the wind speeds 
relative to the storm motion which are the values punched on cards and stored 
on magnetic tape. To convert to wind speed relative to the ground the frontal 
motion must be added to W . The next subsections describe tables and 
graphical illustrations of the wind speed profiles. 

4.1.1 Tables of Wind Speed 

Tables of synoptic cold frontal passage wind speeds are given in 
Appendix 4A. The tables have, due to their length, been split into six parts 
covering grid stations 1 through 20, 21 through 40, etc. The tabulated 
values of W are values relative to the storm motion. The frontal speed is 

A 

given at the top of the table for converting W to its value relative to 

A 

ground. 

The frontal storm case numbers' designation as 1 and 2 in this report 
are listed in the upper left-hand corner. 
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Also, listed at the top of the table is the horizontal length scales for 
the given wind record. The horizontal extent of each data set varies because 
of the data reduction procedure. Hence, the length of field, L, in kilo- 
meters and the horizontal grid spacing. Ax, are specified on each table. The 
vertical extent of each field is taken as 500 m with 50 m vertical grid 
spacing. 

4.1.2 Illustrations 

Illustrations of the longitudinal, lateral, and vertical wind speeds 
encountered along four 3° glide slope through each front are provided in 
Appendix 4B. The glide slopes are displaced in increments of 4.3 km for Case 
1 and 3.3 km for Case 2. Each profile is the wind seen by a airplane travel- 
ing along the flight path drawn across the streamlines as il lustra ted in the 
first figure of the appendix. Note that the streamlines plotted are relative 
to the speed of the front which for reference purposes is indicated with the 
vertical dashed line on the horizontal wind speed profile. The wind speed 
profiles are relative to the fixed earth frame of reference. 

4.1.3 Computer Program 

The card deck or magnetic tape on which the frontal wind speeds are 
tabulated can be directly coupled with the program for thunderstorms given in 
Appendix 2C and requires and provides identical input/output statements. The 
only modification needed is to increase all dimension statements, do-loop 
statements and read-in control integers from 41 to 113. 

The turbulence option can also be called for if desired. The turbulence 
intensities and and length scales Lw^* L^y and Lwz used in the 

program all have values as estimated for thunderstorms, however, in lieu of 
better information on turbulence in major fronts, these values can be used. 

4.2 Turbulence Model 

No data on turbulence in large frontal flows were found in the litera- 
ture. It is therefore suggested that the model proposed for the thunderstorm 
be employed along with the lower values of the input parameters. Recall that 
data on turbulence in thunderstorms is also very limited in the region near 
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the ground, and the thunderstorm turbulence model is itself an approximation. 
Again, it is pointed out that only the high frequency components of the 
turbulence (i.e., frequencies greater than 0.1. Hz) have been lost in the data 
since the quasi -steady wind field given represent 10 sec averages. There is, 
however, a very great need to obtain turbulence data in thunderstorms and 
large fronts in the vicinity of the ground or terminal areas of an airport. 

4. 3 Stati sti cal Model 

Due to the paucity of detailed wind speed data for major fronts, no 
attempt is made to establish a statistical model of the wind fields. 

4.4 Applications of Frontal Wind Shea)^ Model 

The frontal wind shear model is applied exactly the same as the thunder- 
storm models described in Section 2,6. The user, however, has only limited 
choice of wind fields to use in the simulation since only two fields are 
available, Keenan [4-2] has used thunderstorm case numbers 2 and 3 as fronts 
since these two thunderstorms behave in a manner similar to fronts. However, 
they are extremely short in length and are not long enough for a full length 
simulation based on an aircraft approaching from 500 m along a 3® glide 
slope. Such an approach requires a horizontal fetch of 9.5 km. Fortunately 
the frontal data that are available for the two cases of wind shear studied 
extend over a sufficiently long distance that a sequence of one or more 
landings or takeoffs, can be simulated using these data. Up to two landings 
assuming a 3 mile separation distance can be simulated. Moreover, the wind 
fields encountered along flight paths separated by small distances show quite 
dissimilar wind speed profiles as evidenced from the graphical display in 
Appendix 4B. 

Again, it is recommended that turbulence models be used in conjunction 
with the quasi-steady state wind speeds utilized in developing the frontal 
wind shear model. It is anticipated that extreme wind speeds will be somewhat 
averaged out using 10 sec averages, just as described for the thunderstorm 
models, and therefore to achieve a correct simulation high gusts should be 
introduced in the form of turbulence. The turbulence intensities used in the 
model, however, are still open to question, and it is proposed that the lower 
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extremes of turbulent intensities quoted for the thunderstorm cases be util- 
ized for frontal wind shear. 

Much more data on fronts is needed before an adequate model of a synoptic 
front wind field is possible. Such data may be available at the National 
Severe Storms Laboratory in Norman, Oklahoma, but require further processing. 
Also, to utilize the 10 sec average wind speeds generated by the NSSL Program 
effectively, more must be determined about the turbulence levels in major 
fronts . 

One additional set of data that the engineer interested in carrying out 
warm front simulation should be aware of is a one-dimensional profile compiled 
by Keenan [4-2]. This model was developed from data obtained in an accident 
which occurred in Tokyo in 1966. Table 4-1 illustrates the warm front 
one-dimensional model. 


TABLE 4-1 

WARM FRONT WIND SPEED PROFILE [4-2]; W^ = 0 


z (m) 

W^ (m s"^) 

Wy (m s"^) 

6.1 

-6.6 

3.9 

30.5 

-3.6 

-2.1 

61.0 

-4.4 

-12.1 

91.4 

0 

-20.6 

121.9 

5.4 

-30.4 

152.4 

5.9 

-33.5 

182.9 

5.8 

- 32 ! 9 

213.3 

5.8 

-32.9 

243.8 

5.7 

-32.4 

274.3 

5.6 

-31 .9 

304.8 

5.6 

-31.4 

457.2 

5.9 

-31 .4 


86 



4.5 References 


4-1. Goff, R. Craig. "Thunderstorm Outflow Kinematics and Dynamics," NOAA 
Tech Memo ERL NSSL-75, December 1975. 

4-2. Keenan, M. 6. Personal communications, Stanford Research Institute, 
Manlo Park, California, October 1977. 


87 



SECTION 5.0 


SUMMARY 

Mathematical models of wind fields associated with thunderstorms, 
neutral and stable boundary layers, and storm fronts have been developed 
based on reported experimental data. These data have been tabulated and 
presented in Appendices 2A, 3A and 4A, respectively. Computer programs have 
been developed which utilize table lookup and interpolation routines to com- 
pute wind speeds at spatial positions and under stability conditions as 
called for by the main calling computer program. These computer programs are 
presented in Appendices 2C and 3C. 

Turbulence simulation techniques have been developed which permit the 
super-positioning of random fluctuations having the characteristics of 
turbulence associated with the atmospheric phenomena of interest. 

The wind speed models for thunderstorms and major fronts represent the 
three components of wind velocity in a vertical plane 500 m in extent and of 
variable lengths in the horizontal direction. In simulating an approach, 
application of these models assumes that the flight path lies in the plane of 
the wind speeds with that plane being centered over the runway. These models, 
therefore, give realistic simulation of wind shear due to a thunderstorm 
moving parallel to the runway and allow the flight mechanics computation to 
have departures from the flight path in the vertical direction. Any 
departure of the aircraft out of the plane in the lateral direction is not 
accurately modelled. However, evidence is shown that the correlation between 
the wind components in the x- and z-direction is large and, therefore, it is 
anticipated that the same is true in the y-direction. Therefore the wind 
disturbances encountered by an aircraft displaced by the wind a reasonable 
distance out of the plane of the wind field will be nearly similar to the in- 
plane winds. 

To simulate a flight path which is diagonal to the direction of the 
motion of the front, one must assume that the wind fields modelled are two- 
dimensional over a very large lateral scale. This, however, is not likely 
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and to provide wind fields which permit this type of simulation will require 
considerably more research into the three-dimensional nature of thunderstorms 
and major fronts. 

The turbulence model developed to accompany the mean wind field data are 
believed to be quite realistic with the exception that the value of the 
turbulence parameters in the frontal case are vertically unknown. The turbu- 
lence model for the neutral and stable boundary layers is based on extensive 
data sets and theoretical analysis and is believed to be very representative 
of what happens in the actual atmosphere. The thunderstorm turbulence model 
is based on advanced turbulence concepts and does include large-scale velocity 
fluctuations of a non-Gaussian nature, however, coherence matching between 
large gusts which are in the development stage [5-1] should possibly be added 
to the simulation technique. Finally, the parametric inputs to the model 
require better definition which can only be achieved with more experimental 
effort. Currently, it is proposed to use the values of turbulent intensity 
measured at elevations of several kilometers because these are the better 
defined values. Moreover, comparing the thunderstorm data measured by ground 
based towers to those backed out of on-flight data recorders shows that there 
is quite a large difference between 10 sec averaged data and much shorter 
averaging period data. This suggests the available thunderstorm data be re- 
analyzed using a smaller time increment. 

There is considerable need to measure turbulence in the vicinity of the 
ground, particularly in the downdraft center of thunderstorms and also to 
measure the gust gradient which could occur in these severe storms phenomena. 
The same conclusion holds true for large-scale frontal motion where very 
limited, if any, wind velocity measurements as well as turbulence measure- 
ments have been made. The type of measurement which appears most needed 
would be from an array of towers, preferably 500 m in height, and distributed 
in a line at least over 1 km in extent. Such an array of towers is obviously 
very expensive and the data reduction associated with the vast amounts of 
data generated is, in turn, costly. A tower array is available at NASA 
Marshall Space Flight Center with the limitations that the towers are only 
24 m in height [5-2]. 
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As techniques for remote sensing of wind speed become available, typical 
examples being laser Doppler, acoustic Doppler or Doppler radar methods, more 
probing of thunderstorms and of major synoptic fronts will undoubtedly 
increase. 

Aircraft are also being used to study thunderstorms and are providing 
meaningful data at great altitudes. It is not likely that much probing of 
thunderstorms or severe frontal weather will be carried out with aircraft 
below 500 m. However, data at this level is extremely critical in developing 
simulator wind field models for flight/hazard definition. 

A possible technique for obtaining the needed data would be to couple 
aircraft measurements with those of the tower array currently in existence at 
the Atmospheric Sciences Division, NASA/Marshall Space Flight Center, 
Huntsville, Alabama. A coordinated program utilizing aircraft passes over 
the tower array with the towers in operation would map out a well defined 
two-dimensional wind field and could also provide three-dimensional wind 
fields when the angle of the approaching storm is oblique to the tower array. 

Mathematical models currently being developed have promise for use in 
wind shear simulation applications, however, they are, in turn, dependent 
upon input parameters which must be experimentally measured before meaningful 
results can be calculated. 

A statistical model for the stable and neutral boundary layer has been 
developed which allows the probability of a given wind speed and stability 
condition occurring simultaneously to be estimated. The model also predicts 
the frequency of the combined conditions occurring within a specified period 
of time. This model assumes a Rayleigh distribution of wind speeds for any 
given site in the mainland U.S.A. This is a justifiable assumption based on 
a study of 138 sites by Justus, et al . [5-3]. The probability distribution 
of Richardson numbers is not quite as well validated, being based on only two 
sets of data. It is believed, however, that this is a good estimate of the 
distribution of stability conditions with wind speeds. 

The statistical model for the thunderstorm is rather weak and for the 
synoptic fronts nonexistent. The best estimate of the probability of a given 
magnitude of wind shear in a thunderstorm that is currently available is 
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based on the ensemble average of 20 thunderstorm cases along a prescribed 
flight path and the standard deviation from the averaged value. These com- 
puted standard deviations are added at each point along the flight path as if 
they were independent points. This is a questionable assumption and further 
work is needed to develop the statistical model of a thunderstorm. 

Before a statistical model of fronts can be developed much more data is 
needed. Due to the lack of measurements of detailed wind speed profiles 
through large fronts and also, to an extent, through thunderstorms, reliable 
simulation flight/hazard wind speed models cannot be developed. 

The conclusion of this report is that mathematical models of wind shear 
in thunderstorms, stable and neutral boundary layers and fronts for flight/ 
hazard simulation studies have been developed which, based on state of the 
art knowledge, allow realistic simulation of flight through hazardous 
atmospheric wind shear. These models, however, represent a very good 
beginning but they should be continuously updated as information becomes 
available and standardized to establish consistent criteria for pilot 
training, avionics development and aircraft design and certification. 
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APPENDICES 

The following appendices are numbered according to the section of the 
main body of the report to which they pertain. For example. Appendix 2 
followed by an alphabetical symbol contains information relating to Section 2, 
Thunderstorm Wind Shear. Appendices 3 and 4 are referenced in the same 
format. 
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LIST OF ABBREVIATIONS AND SYMBOLS 

Friction velocity 

Stability parameter, KU*/fL' 

Coriolis parameter 

Monin-Obukhov stability length scale 

Horizontal extent of thunderstorm data set 

von Karman's constant taken as 0.4 

Length scale of longitudinal velocity fluctuations 

Length scale of lateral velocity fluctuations 

Length scale of vertical velocity fluctuations 

Mean speed of storm front 

Longitudinal wind speed 

Lateral wind speed 

Vertical wind speed 

Fluctuating longitudinal wind speed 

Fluctuating lateral wind speed 

Fluctuating vertical wind speed 

Surface roughness 

Dimensional height 

Dimensionless height z = zf/u* 

Rossby number u*/fz^ 

Power density spectrum function 
Cyclic frequency 

Turbulence intensity of a wind speed fluctuation 
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Statistical correlation between a and 0 wind speed components 

Ensemble average 

Denotes dimensionless parameter 



APPENDIX 2A 


TABULATED THUNDERSTORM DATA 

Appendix 2A contains the tabulated data interpolated from the thunder- 
storm wind speed data presented and documented in Goff [2-1]. These data are 
the fundamental data sets upon which the wind shear model for thunderstorms 
is based. The user can compute wind speed profiles either by hand directly 
from the tables or by coupling the data with the lookup computer program 
provided in Appendix 2C. 

The rows and columns of the table are numbered from the left-hand bottom 
corner. The numbers increase from left to right in the positive x-direction 
of the original data and from bottom to top in the positive z-direction of 
the original data. The wind speeds are given in units of meters per sec, 
m s'\ with W being positive in the direction of frontal motion, Vi being 

A Cm 

positive upward and being positive into the plane of the paper. The 
upper portion of the table covers columns 1 through 21 and the lower table 
covers columns 21 through 41. (Note column 21 is repeated for synmetry and 
clarity of presentation.) 

The thunderstorm case numbers' designation for this report are listed 
at the top of the table. The letters in parentheses and the following series 
number correspond to the thunderstorm designation given by Goff [2-1]. 

Also listed at the top of the table are the frontal speed, W , and the 

A 

horizontal length scales for the given wind record. The horizontal extent 
of each data set varies because of the data reduction procedure. Hence, the 
length of field, L, in meters and the horizontal grid spacing. Ax, (m) are 
specified on each table. The vertical extent of each field is taken as 
500 m with 50 m vertical grid spacing. Figure 2A-1 illustrates the identi- 
fication format used for all tables. 
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TABLE DESIGNATING SYMBOLS 
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Table 2A-1(a) Wind Speed in Direction of Frontal Motion, 

Case 1(A) Series No. 0446, 

6*1 m s \ "^ength of Field L = 4000 m/ Horizontal Increment Ax = 100 meters 
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Table 2A-l(b) Wind Speed Perpendicular to Frontal Motion, 

Case 1(A) Series No. 0446, 

1.1 m s Length of Field L = 4000 m Horizontal Increment Ax = 100 meters 


vO 

U 
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Table 2A-l{c) Vertical Wind Speed, W^ 

Case 1(A) Series No. 0446, 

Length of Field L = 4000 m/ Horizontal Increment Ax = 100 meters 
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Table 2A-2(a) Wind Speed in Direction of Frontal Motion, W 

X 

Case 2(B) Series No. 1314, 

0ms , Length of Field L = 3278,8 m. Horizontal Increment Ax = 81.97 meters^ 
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Table 2A-2(c) Vertical Wind Speed, 

Case 2(B) Series No. 1314, 

0 n\ s , Length of Field L = 3278.8 m. Horizontal Increment Ax = 81.97 meters 
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Table 2A-3(a) Wind Speed in Direction of Frontal Motion, 

Case 3(C) Series No. 1731, 

Length of Field L = 5714.4' m# Horizontal Increment Ax 142.86 meters 
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Table 2A-3(b) Wind Speed Perpendicular to Frontal Motion, 

Case 3(C) Series No. 1731, 

8.6 m s~^. Length of Field L = 5714.4 ra, Horizontal Increment Ax « 142.86 meters 
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Table 2A-3(c) Vertical Wind Speed, 

Case 3(C) Series No. 1731, 

8.6 m s~^. Length of Field L = 5714.4. m. Horizontal Increment Ax = 142.86 meters 
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Table 2A-4(a) Wind Speed in Direction of Frontal Motion. W 

Case 4(D) Series No. 1459, 

6 m s"V Length of Field L = 7692.4 m. Horizontal Increment 
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Table 2A-4(b) Wind Speed Perpendicular to Frontal Motion, 

Case 4(D) Series No. 1459, 

11.6 m s , Length of Field L = 7692.4 m, Horizontal Increment Ax = 192.31 meters 
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Table 2A-4(c) Vertical Wind Speed, 

Case 4(D) Series No. 1459, 

11.6m s , Length of Field L = 7692.4 m. Horizontal Increment Ax = 192.31 meters 
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Table 2A-5(a) Wind Speed in Direction of Frontal Motion, 

^ Case 5{E) Series No. 1924, 

=16.7 m s , Length of Field L = 11428.4 m. Horizontal Increment Ax » 285.71 meters 
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Table 2A“5(b) Wind Speed Perpendicular to Frontal Motion, 

Case 5(E) Series No. 1924, 

16.7m s , Length of Field L = 11428.4 m , Horizontal Increment Ax = 285.71 meters 
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Table 2A-5(c) Vertical Wind Speed, 

Case 5(E) Series No. 1924, 

16.7 m s~^. Length of Field L * 11428.4 Hi, Horizontal Increment Ax =» 285.71 meters 
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Table 2A-6(a) Wind Speed in Direction of Frontal Motion, 

Case 6 (f) Series No. 1933, 

11.0 m s ^ Length of Field L = 7272.8 m. Horizontal Increment Ax = 181.82 meters 
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Table 2A-6(b) Wind Speed Perpendicular to Frontal Motion, W 

Case 6(F) Series No. 1933, 

11.0m s , Length of Field L = 7272.8 m. Horizontal Increment Ax = 181.82 meters 
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Table 2A-6(c) Vertical Wind Speed, 

^ Case 6(F) Series No. 1933, 

11.0 m s , Length of Field L = 7272.8 m. Horizontal Increment Ax = 181.82 meters 
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Table 2A-7(a) Wind Speed in Direction of Frontal Motion, 

Case 7(G) Series No. 1942, 

11.8m s~^ Length of Field L = 7843.2 m. Horizontal Increment Ax =196.08 meters 





00 

00 

oo 

CD 

9 

«4 

9 

9 

9 

9 

d 



>«■ 

d 



d 

d 

d 

d 

d 

9 


» 

, 

• 

• 

• 

.« 

• . 

»■ 

« 

• 

• 


• 

• 

• 


• 

• 

• 

• 

• 

ft 

ft 


o» 

o 


*14 


o 

9 

lA 

d 

O 

9 

-<• 

9 

9 

9 

A 

A 


d 

d 

. 44 

O 

O 


<r* 

d 

d 

d 

d 

d 

44 

4>t 

44 

44 



• 

1 

1 

1 

a 

a 

a 

a 

a 

a 



■■ 

♦ 

o 

m 

9 

d 

9 

9 

9 

A 

9 


d 

9 

r«- 

9 

9 

4T 

o 

A 

9 

44 

d 

o 

• 


• 

• 


.• 

• 

• 

• 

• 


o 

• 

• 

• 

• 

• 


• 


• 

m 

ft 

Of 

O 


d 

d 


o 

9 

un 

d 


9 


9 

9 

r- 

A 


d 

d 


O 

o 

44 



d 

d 

d 

d 

d 

■•4 

«4 

44 

44 



1 

1 

t 

a 

a 

a 

a 

a 

a 





o 

d 

9 

9 

9 

9 

CD 

o 

d 

9 


A 

d 

d 

d 

d 

.d 

9 

9 

A 

A 

9 

Ov 

• 


• 

* 

• 

• 

*. 

• 

• 

• 

.• 

9 

• 

■.• 

• 

• 

• 


• 

• 

• 

A 

ft 



d 

d 

d 

44 

O 

9 

lA 

A 

d 

9 

d 

9 

9 

.d 

A 


d 

O 

O 

O 

o 

44 



d 

d 

d 

d 

d 


.44 

44 

44 



f 

• 

1 

a 

a 

a 

a 

a 

a 




m 

4» 


m 

lA 

9 

9 

9 

d 

44 

9 


.44 

A 

.p= 

A 

9 

9 

d 

d 

d 

9 

9 

00 


• 


m 

• 

• 

.• . 


• 


• 

9 

• 

« 

• 

• 

. • 

• 

* 

* 

» 

ft 

ft 


d 

d 

d 

d 

•4 

O 

9 

aO 

9 

d 

O 

d 

9 

d 

A 


d 

44 

o 

9 

o 

O 

d 


d 

d 

d 

d 

d 

d 



44 


94 


• 

1 

a 

» 

a 

a 


a 

a 




♦ 

d 


9 

o 

44 

9 

9 


9 

9 


A 

9 

d 

.O 

d 

d 

A 

A 

d 


9 

»««• 

• 

• 

• 

.• 

• 

• 

• 

• 

• 

• 

.« ' 

d 

« 

• 

• 

.• 

• 

• 

* 

• 

ft 

.« 

ft 




d 


«k4 

«4 

9 

d 

d 

d 

d 

d 

r- 

A 

A 

A 

d 


44 

■o 

44 

d 

d 


d 

d 

d 

d 

d 

d 

44 

44 

44 

44.. 

•4 


1 

9 

• 

a 

1 

a 








00 

oo 

9 

d 

d 

d 

9 

9 

9 

■X 


9 

d 

.O 

A 

■P4. 

9 

A 

.44 

9 

d 

9 


m 

• 

m 

• 

• 

• 

.• 

• 

• 

• 


A 

• 

• 

« 

• 

9 

. 9 

« 


.9 

ft 

ft 



“O 

d 

d 

.•4 

«4 

o 

9 

9 

A 

A 

d 

A 

A 

A 

d 

d 

^ o 

d 

44 

d 

d 

A 


d 

d 

d 

d 

d 

d 

d 


.. .44 

44 

«4 


1 

I 

• 

1 

a 

• 








tn 

o 

PA 

d 

9 

d 

9 

9 

44 

9 


A 

.44 

9 

9 

d 

d 

9 

» 

d 

A 

9 


• 

• 

• 

. • 


• 

• 

• 

9 

• 

• 

A 

• 

• 

« 

• 

• 

9 

« 


• 

ft 

ft 

«-4 


-n 

d 

d 


*4 

.O 

9 

9 

CD 

9 

d 

A 

A 

d 

d 

d 

o 

d 

44 


9 

A 


<>• 

d 

d 

d 

d 

d 

d 

d 


44 , 

<4 


1 

1 

1 

1 

a 










«<■ 

tA 

o 

9 


9 

■4' 

d 

9 


O 

9 

A 

d 

d 


o 

A 

9 

9 

9 


• 

• 

• 

• 

• 


• 

4 

.9 


« 


• 

• 

• 

« 

m 

• 

• 

• 

9 

• 


.«<i4 

m 

rn 

d 

d 

d 

•4 


O 

.O 

o 

9 

d 

A 

d 

d 

d 

o 


A 

d 

A 

A 

A 


d 

d 

d 

d 

d ^ 


d 

d 

d 

d 



1 

• 

a 

a 

a 









00 

d 

o 

9 

d 

jn 

d 

9 

44 

9 


d 

00 


d 

9 

9 


94 

d 

.44 

9 

m 


• 

• 

• 


M 

• 

• 


• 

• 

d 

• 

• 


.» 

• 

> 

9 



• 

• 

•t 


d 

d 

d 

•4 

.44 

o 

o 

9 

9 

9 

d 

d 

d 

«4 

o 

O 

o4 

A 

A 

d 

d 

f-x 


d 

d 

d 

d 

d 

d 

d 

d 

4i| 

•4 ■ 

*4 


1 

» 

• 

a 










lA 

o 

4T 


9 

9 

9 

d 

9 

.9 


d 

A 

d 

9 

■ d 

o 

d 

9 

CO 

d 

9 

:>! 


• 

• 

• 

• 


V 

• 

.» 

« 

.• 

d 

m 


• 


■m 

« 


• 

9 




d 

d 

d 

•4 

«4 

O 

9 

9 

9 

d 

d 

d 

d 

44 

o 

O 

d 


r» 

A 

9 

9 

9 


d 

d 

d 

d 

d 

d 

«4 . 

i 44 

44 


.44 


1 

• 

1 










d 

«<« 

00 

9 

d 

d 


9 

d 

A 

9 


9 

A 

44 

d 

9 

d 

A 

A 

d 

A 

9 


• 

. • 


• 

.• 

» 

• 

• 

' • 

• 


*4 


'• 

m 

9 

• 

• 


9 


• 

ft 


d 

d 


O 

o 

O 

9 

9 

d 

d 

A 

d 

O 

O 

•4 

d 

d 

A 

9 

9 

o 

9 

9 


d 

d 

d 

d 

d 

d 

•4 

44 

«4 

•4 

•4 


1 

1 







44 




m 

00 

9 

d 

d 

9 

o 

.o 

4T 

O 

A 


.A 

9 

A 

9 

d 


A 

44 

9 

CP 

9 

o 


• 

• 

*■ 

e 

• 

• 

• 

.• 

• 


o 

• 

• 

« 

• 

« 

• 

• 

.• 

ft 

ft 

ft 

•-« 



o 

o 

O 

9 

9 

9 

d 

d 

A 

d 

o 

O 

d 

d 

A 

9 

9 

o 

44 

o 

O 


d 

d 

d 

d 

d 

«4 

44 

.44 

44 

■44 

»4 









.44 

44 

d 

44 




9 

CD 

30 

*4 

9 

vO 

d 

A 

d 


9 

9 

9 

9 

9 

A 

9 

9 

d 

A 

9 


• 

• 

* 

• 

• 

• 

• 

• 

« 

• 


9 

« 

• 

• 


• 

• 

• 

• 

.ft 

• 

ft 



o 

o» 

9 

9 

9 

9 

9 

d 

A 

A 

d 

44 

d 

d 

A 

A 

O 

O 

44 

d 

d 

.44 


d 

d 



.■4 

44 

44 

44 

44 

44 

44 







44 

■44 

44 

d 

9 

44 


<n 

00 



CD 

(A 

d 

d 

9 

d 

9 


9 

9 

9 

9 

9 

9 

44 

9 

9 

d 

9 

3> 

• 

• 

• 

• 

• 

• 


.• 

• 

• 

• 

9 

• 

• 

• 

• 

• 

• 

• 

.• 

ft 

ft. 

.* 


o 

a* 


9 

9 

9 

9 

9 

A 

A 


d 

A 

A 

'O 


9 

d 

d 

d 

A 

A 

d 


d. 

•4 


•■4 

4« 

.«4 

44 

.44 

■ 44 

•4 

.44 






44 

44 

44 

44 

•4 

44 

44 




vA 


00 

9 

9 

9 

9 

9 

.d 


9 

9 

9 

9 

9 

44 


9 

d 

9 

9 

r- 

> 


• 


* 


• 


• 

*' 

:» 

d 

• 

• 

• 

.• 

.• 

* 



.ft 

ft 

ft 



O' 



d 

d 

d 

d 

A 

A 

4** 

d 

r*. 

9 

9 

d 


4T 

A 

A 

A 

A 

A 



«<4 

«>4 

-* 

•4 


*4 

44 

•4 

.«4 

.44 





.44 

44 

44 

44 

.44 

44 

44 

44 


00 

vO 

00 

9 

9 

9 

sA 

*n 

d 

A 



9 

9 

9 

o 

»>• 

A 


9 

CO 

A 

9 

^o 


« 

• 

• 

• 

• 



* 

• 

« 

A 

• 

• 

.« 

• 

.• 

• 

• 

■• 

ft 


ft 


00 

00 



d 

d 


d 

A 

A 


d 

O 

d 

d 

'<r 

A 

A 


d 

f*> 

d 

A 



•M. 


v4 

44 

.-4 

44 


■.-4 

44 

44 


44 

44 

44 

•4 

44 

44 

.44 

44 

44 

.44 

44 


CD 


CD 

9 

in 

/4 

d 

4-» 





9 


9 ' 

•4 

A 

9 

d 

d 

o 

d 

9 

(A 

». 

• 


• 

• 

* 

• 

• . 

.• . 

4 

• 

m 

• 

« 

• 

• 

• 

• 

»■ 

.« 

ft 


ft 


■r^ 

P»- 


.r» 

d 

d 

d 

d 

A 

A 


d 

d 


A 

n 

A 

A 

xO 


A 

A 

A 




-4 

•4 

44 

•4 

•4 

44 

44 

•4 

— 


•4 

44 

-.44 

.44 

44 

44 


44 

44 

44 

44 


Gh 

d 

CD 

r- 


•4 

9 

9 

9 

9 

CD 


9 

9 

X 

•4^ 


X 

d 

A 

,.4 

9 

9 



• 

• 



• 


• 

• 



-*• 




• 


• 



ft 




.1^ 

r- 


T- 

d 

d 

9 

Vi' 

A 

4T 

d 

CN 


A 

A 

A 

A 


r- 

A 

4f 

d 

d 





•4 

•4 

•4 


. 

44 

44 

44 


44 




4« 

44 

.44 

-4 

44 

44 

.44 




CD 


d 

no 

9 

A 

A 

A 

A 


9 

A 

A 

9 

9 

9 

44 


o 

.44 

9 


• 

• 


.« 

» 

• 



• 

• 

.• 

d 

• 

« 

• 


.• 

• 

■ 

• 

ft 

ft 



x> 


,r«- 

r- 

:d 

x> 

9 

9 

A 


d 

d 

A 

'4* 

ro 

A 

A 

9 

OP 

•A 


d 

9 


*>4 




«4 

•4 

44 


»4 

44 

44 



44 


■-^ 



.44 

44 

-4 

44 



d 



Xt 


A 

d 

•4 

•4 


.44 


JO 

d 

»4 

9 

9 

9 

o 

» 

o 

A 

d 


• 




• 

• 

». 


.• 

• 

* 

d 

;• 

• 

• 

• 


.» 


ft 

ft 

ft 

ft 


CD 

r* 




A 

-6 

A 

A 


d 

d 

*>4 

9 

o 

9 

9 

9 

9 

A 

9 

ft4 

9 





-4 

•4 

44 

.4-4 

44 

-.44 

•4 

44 




d 

44 


•4 

.44 


44 

44 



•*» 

•n 

o 

OD 

9 

d 

o 

CD 

9 

CO 

9 


jp 

9 

9 

9 

9 


9 

9 

9 

A 

d 


• 

• 

• 



« 


:• 

• 


• 

44 



* 

• 

• 



ft 


ft 

ft 


•c 

T 

•» 

r* 

9 

9 

A 

A 

rr 

d 

d 

d 

9 

O 

•4 

.44 

44 

o 

9 

•A 

d 

O 

9 





•4 

■4 

4«4 


.4.1 

•4 

44 




d 

d 

d 

d 

d 

44 

.44 

44 

44 




O 

9 

•o 


9 

lA 

A** 

d 

d 

•4 


44 

o 

9 

9 

4. 

A 

A 


d 

d 




115 



c 

o 


■M 

O 


<C 

4 -> 

d 

o 

s- 




u 

•f— 

XJ 

c 

QJ 

O- 

S-. 

OJ 

O- 

-o 

0 ) 

a. 

CO 

XJ 

c 


CO 

M 

I 

00 

o 

vO 

<T\ 


X 

o 

+J 

a 

<u 

o 

c 


<d 

-p 

G 

o 

U 

^ o 

fN JU 

iH a 




<U 

CO 

<T 3 

U 


rsi 


O 

cn 

m ^ 
0) 00 
•H 

u 

<U II 


X 

f— i 

Q) 


I 

<c 

CNJ 

<u 

r— “ 

XI 

ftJ 


x: 

tn 

G 

(U 

h:i 


CO 

a 

00 


rH 



116 



Table 2A-7{c) Vertical Wind Speed, 

^ Case 7(G) Series No. 1942, 

11.8 m s , Length of Field L = 7843.2 m , Horizontal Increment Ax=l96o08 meters 


ORIGINAL PAGE B 
OF POOR QUALrra 


II 


Is 


X 



117 



3: 


A 


c 


o 


•r~ 


4-> 


O 


2: 




fO 


-M 


C 


O 


$- 


Ll- 

m 

CN 

if- 

•H 

O 




c 


o 

• 

•r— 

O 

4-> 


O 


O) 

CO 

i- 

0) 

*r— 

•H 

Q 

f-j 


(U 

C 

CO 

•1- 



N 

<U 

rc 

0) 


CL 

00 

tn 

cu 

“O 

CO 

c 

d 

•r~ 

o 

3K 



I 

<C 

CSJ 

<u 


fO 


CO 

<D 

(D 

6 

as 

00 

00 

m 


<1 


c 

(U 

e 

OJ 

u 

o 

d 


d 

o 

N 

•H 

U 

o 


so 

in 

m 

m 

li 

T3 

rH 

Qj 

•H 

4-1 

O 

4J 

W) 

C 

0) 


w 

a 

m 

00 




118 



Table 2A-8(b) Wind Speed Perpendicular to Frontal Motion, 

Case 8(H) Series No. 1712, 

=8.5 m s~^ Length of Field L = 5555.6 m , Horizontal Increment Ax = 196.08 meters 
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Table 2A-8(c) Vertical Wind Speed, 

Case 8(H) Series No. 1712, 

Length of Field L = 5555.6 m. Horizontal Increment Ax = 138.89 meters 
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Table 2A-9(a) Wind Speed in Direction of Frontal Motion, 

Case 9(1) Series No. 1507, 

11.5 ms, Length of Field L = 7407.6 m. Horizontal Increment Ax = 185.19 meters 
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Table 2A-10(a) Wind Speed in Direction of Frontal Motion, 

Case 10 (J) Series No. 2206, 

13.1m s~^ Length of Field L = 8510.8 m.. Horizontal Increment Ax = 212.77 meters 
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Table 2A~10(b) Wind Speed Perpendicular to Frontal Motion, 

Case 10 (J) Series No. 2206, 

13.1 m s ^ Length of Field L = 8510.8 m, Horizontal Increment Ax = 212.77 meters 
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Table 2A-10(c) Vertical Wind Speed, 

Case 10 (J) Series No. 2206, 

13.1m s~^. Length of Field L = 8510.8 m. Horizontal Increment Ax = 212.77 meters 
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Table 2A-n(a) Wind Speed in Direction of Frontal Motion, 

Case 11 (K) Series No. 2118, 

12.4 m s~^ Length of Field L = 8163.2 ra , Horizontal Increment Ax = 204.08 meters 


ORIGINAL PAGl^ 
OF POOR QUAILH^i 


fl 

1:5 
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Table 2A-n(b) Wind Speed Perpendicular to Frontal Motion, 

Case 11 (K) Series No. 2118, 

12.4 m s~\ Length of Field L = 8163.2 m , Horizontal Increment Ax = 204.08 meters 




l3s 
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Table 2A-ll(c) Vertical Wind Speed, 

Case IKK) Series No. 2118, 

12.4 m s~^» Length of Field L = 8163.2 Horizontal Increment Ax = 204 .OS' meters 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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Table 2A-12(a) Wind Speed in Direction of Frontal Motion, 

Case 12 (L) Series No. 1759, 

m s"^ Length of Field L = 3636.0 m , Horizontal Increment Ax = 90.90 meters 



130 


7.S 7.8 6.7 5.5 i.4 i.O 2.6 2.3 1.9 1.5 0.8 0.0 -O.? 



ORIGINAL PAGE IS 
OF POOR QUALITY 

W 



0 

+J 



' 

: 






• 


- — " • 














0 


m 

o 


O 

o 

CM 

^ftlft 

<M| 

•M 

•ft 

O 


«i4 

•ft 

•ft 

ftftt 

Ok 

M* 

9 

•ft 

•ft 


e 


6 

it 

• 

ft 

,ft 

o 

ft 

o 

o 

ft' 

o 

.ft 

o 

ft 

o 

ft 

o 

ft 

o 

ft 

9 

HT 

ft 

•ft 

ft 

•ft 

« 

• 

o 

• 

o 

-• 

o 

* 

o 

ft 

9 

• 

o 

• 

o 









a 

a 

a 

a 

a 










a 

a 

a 



O 


























c^ 


o 

o 


na 

CM 

CM 

CM 


'•M 

•ft 

-9 


© 

o 

o 

o 

Ok 

*o 

O 

•ft 

•ft 

•ft 

o 


o 

• 

ft 

ft 

• 

ft 

ft 

ft 

ft 

ft 

• 

ft 

9 

# 

• 

• 

ft 

• 

• 

ft 

• 

ft 

• 

• 


• 

<N 

o 

o 

o 

o 

o 

O 

O 

o 

O 

o 

O 

HT 

•ft 

•ft 


•ft 

9 

9 

9 

9 

o 

9 

9 


o 




1 

a 

a 

a 

a 

a 

a 

a 










1 

a 

a 



<T> 


w 

fM 

ro 

m 

cn 

m 

CM 

•f* 


.•ft 

o 


o 

o 

o 

o 

9 

•O 

9 

9 


•ft 

O 



0^ 

■m- 

-ft 

ft 

' • 

ft 

. -ft 

ft 

ft 

ft 

• 

.ft 

Ok 

• 

ft 

• 

• 

* 

ft 


ft 

• 

• 

• 


il 

«•« 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

9 

m 

•ft 

•ft 

•ft 

•4 

9 

9 

O 

9 

9 

o 

9 

m 



f 

• 

a 

• 

a 

a 

1 

a 

1 

a 











a 

1 


c 

X 


m 


X) 

in 


M* 

m 

CM 

CM 

CM 

9 


Ok 

Ok 

Ok 

o* 

CD 

m 



9 


O 

o 


OD 

-• 

ft 

ft 

ft 

. ft 

ft 

-ft 

ft 

ft 

-• 

- •- 

ao 

• 

• 

• 

» 

ft 

• 

» 

ft 

• 

ft 

• 

•r— 



o 

o 

o 

o 

o 

O 

o 

© 

O 

9 

9 


9 

9 

9 

o 

O 

9 

9 

9 

9 

9 

o 




• 

• 

1 

a 

a 

a 

a 

• 

a 

-a 












1 


o 

-M 



r*' 

r» 


in 

ftf 

m 

CM 

CM 

CM 

9 


Ok 

C9k 

Ok 

OI 

00 

in 

CM 

•ft - 

9 

O 

O’ 

s: 

C! 

1-“ 

• 

• 

ft 

ft 

ft 

» 

ft 

ft 

ft 

ft 

ft 


• 

• 

ft 

• 

ft 

• 


ft 

ft 

• 

ft 

0 


o 

O 

o 

o 

o 

o 

o 

O 

9 

O 

O 

f*3 

O 

O 

9 

o 

o 

9 

O 

O 

9 

9 

9 


0 


• 

• 

1 

t 

a 

a 

a 

t 

1 

a 














<o 

0 


o 

o 

9k 

3> 


«n 

Hi* 

CM 

CM 

CM 

9 



•ft 

9 

o 

9 

kO 

kO 


9 

© 

o 



• 

ft 

ft 

ft 

ft 

-ft 


ft 

• 

-• 

ft 

9 

ft 

• 

ft 

ft 

# 

ft 

ft 

• 

• 

• 

• 

4-> 

w 

•i* 


•M 

O 

O 

o 

' o 

o 

O 

O 

9 

O 

fn 

•ft 

•ft 

•ft 

•ft 

9 

O 

9 

o 

9 

9 

o 

C 

o 


• 

• 

1 

a 

a 

a 

1 

• 

a 

a 














o 

0 




*■< 

o 

3B 

.Ml 

H*» 

CM 

•M 

•ft 

-9 


CM 


•ft 

Ok 

30 

9 

tn 

o 

9 

9 

o 

S- 

Mf 

vn 

• 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

'ft 

• 

m 

ft 

• 

ft 

• 

:ft 

• 

• 

• 

• 

-• 

ft 

U- 


«•* 

«•« 

•4 

w 

•H 

O 

© 

o 

O 

9 

9 

O 

m 

•ft 

•ft 

•ft 

o 

9 

9 

9 

o 

O 

9 

9 


iH 


1 

I 

a 

a 

a 

a 

a 

a 

a 

a 














o fJJ 

0 


fN 


o 

Ok 

kO 

m 

cn 



•ft 

9 


ftf 

M* 

O 

00 

Ok 

r* 

9 

o 

o 

9 

9 

to 

+J 


ft 

ft 

ft 


ft 

.ft 

ft 

ft 

ftl 

.* 

• 

M* 

• 


• 

• 

-• 

-• 

* 

ft 

• 

ft 

• 

J- 

G 

0 


t 

1 

a 

o 

a 

o 

a 

o 

a 

o 

a 

9 

t 

9 

a 

o 

a 

O 


•ft 

•ft 


9 

9 

9 

9 

o 

9 

9 

9 

. 

N 



ftft 

<^ 


»n 


CM 


«M 

^•ft 

9 



in 

9 

Ok 

Ok 

CD 

1- 

•ft 

•ft 

O 

O 

*H 


ft 

ft 

ft 

ft 

ft 

ft 

•' 

ft 

ft 

-• 

. ft 

»n 

• 

• 

• 

• 

• 

• 

ft 

ft 

ft 

• 

ft 

3 S 

in 


«*« 

1 

1 

o 

a 

o 

a 

o 

a 

O 

a 

o 

a 

O 

a 

o 

a 

o 

a 

9 

cn 

•ft 

•ft 

•ft 

o 

9 

9 

O 

9 

O 

o 

o 

o 2 

o 

ffi 



o 


kO 

m 

CM 

CM 

• ^ 


•ft 

9 


00 

r** 

CM 

o 

9 

Ch 

Ok 

m 

•ft 

•4 

9 


fM 

ft 

• 


• 


ft 

ft 

-• 

ft 

• 

ft 

CM 

• 

-• 

ft 

• 

• 

« 

• 

ft 

• 

• 

ft 

T? ? 





o 

o 

o 

O 

O 

o 

9 

o 

O 


•ft 

•ft 

•ft 


•ft 

o 

O 

9 

9 

9 

O 

£= 2 
<u n 

€ 


• 

a 

o 

a 

sO 

a 

tn 

a 

CM 

a 

*•« 

a 

a 

o 

9 

a 

9 

o 


Ok 

Oi 

m 

m 

O 

Ok 

9 

M* 

CM 


o 




ft 

ft 

* 

ft 

ft 

ft 

ft 

ft 

-ft 

.ft 

ft 



t 


• 

• 

• 

• 

ft 

ft 

ft 

m 

sT 


*»« 



o 

o 

O 

o 

o 

o 

9 

9 

9 

m 

•ft 

•ft 

•ft 

•ft 

•ft 

9 

•ft 

O 

O 

9 

9 

5 w 

o 


1 

1 

a 

a 

a 

a 

I 

















Q- 

• 



r>» 

in 


o 

o 

o 

o 

O 

'© 

9 


•ft 

o 

CP 

m 

c*» 

9 

' ftft 

9 

cn 

CM 

O 


3636 

O 

• 

• 

• 

• 

.ft 

ft 

ft 

ft 

ft 

ft 

ft 

9 

ft 


ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

“O S 


*-• 

1 

0 

1 

9k 

0 

1 

o 

a 

M 

o 

o 

o 

«H 

0 0 

o 

«<« 

9 

*<* 

9 

•ft 

O 

9 

cn 

CM 

tn 

CM 

ftf 

•ft 

o 

•ft 

00 

m 

•ft 

m 

•ft 

CM 

9 

o 

M* 

9 

CM 

O 

© 

<U 


• 

• 

• 

• 

ft 

ft 

ft 

ft 

ft 

ft 

'ft 

Ok 

ft 

ft 

ft 


ft 

ft 

ft 

ft 

ft 

• 

ft 

CL 



«H 

O 

o 

O 

o 

o 

o 

9 

O 

O 

9 

CM 

CM 

CM 

CM 

•ft 

•ft 

•ft 

•ft 

9 

9 

o 

p 

^ 0 

I! 


a 

• 

a 

a 






















o 

OD 

r»i 


o 

P4 


CM 

CM 


9 


Ok 


in 

CM 

CD 

m 

M* 

Ok 

M» 

m 

9 

T3 5 


ao 

m 



* 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

9 

ft 

ft 

.ft 

ft 

ft 

ft 

ft 

ft 

ft 

S 

-'• 

c ^ 
U 



a 

0 

1 

o 

a 

o 

a 

o 

© 

o 

9 

9 

9 

© 

CM 

CM 

CM 

CM 

CM 

•ft 


•ft 

9 

9 

9 

'3H 

irH 


o 

•30 


:0 

{M 

•*> 


CM 

•M 

•ft 

9 



•ft 

Ok 

jn 

•ft 

9 ' 

in 

9 

9 

•n 

9 


0 

•H 


• 

m 

» 

• 

ft 

.ft 

ft 

ft 

ft 

-ft 

.ft 

r** 

ft 

ft 

ft 

ft 

• 

ft 

ft 

ft 

ft 


• 



f 

O 

a 

0 

1 

o 

o 

o 

6 

9 

9 

.9 

O 

CM 

fn 


CM 

CM 

CM 


•ft 

•ft 

9 

o 

9 




o 

r*i 

o 

CM 

••• 

m 




9 

9 


CM 

o 


9 

© 

M» 

9 

r- 


CM 

9 



>o 

• 




ft 

ft 

ft 

ft 

ft 

ft 

ft 

9 

ft 


ft 

ft 


ft 

• 

ft 


ft 

ft 




O 

o 

& 

o 

© 

o 

9 

O 

9 

9 

CM 

»n 

f»l 

CM 

CM 

•ft 

•ft 

•ft 

9 

9 

9 

9 

csi 

0 


a 

I 






















r— * 



■■■• 

—» 

<r 

vn 

M3 

© 

CM 

«-« 

9 

•ft 

9 


CM 

o 

O 

m 

•ft 

Ok 

in 

<*»» 

CM 

•ft 

9 

1 



* 

• 

• 

• 

ft 

ft 

ft 

ft 

ft 

-ft 

ft 

in 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

< 

CVJ 

+J 

tJ' 


o 

a 

O 

o 

o 

o 

© 

o 

O 

O 

O 

a 

O 

CM 

m 


CM 

•ft 

•ft 

o 

9 

O 

O 

O 

O 


G 



VD 

00 

r-* 

00 

-r* 

CM 

9 


-•ft 

O 


•• 

o 

O 

CM 

no 

•r 

O 

O 

o 

o 

o 

<D 

0 

hI 


• 

• 


• 

ft 




ft 

ft 

ft 

Hf 

• 


4 


ft 

ft 

ft 

ft 

ft 


ft 


o 

O 

o 

o 

© 

o 

o 

© 

© 

a 

o 

a 

e 

CM 


ftl 

CM 


9 

o 

o 

9 

9 

o 

o 

JCk 



o 


<N 

o 

o 

tc 

Hi* 

CM 

9 

•ft 

9 



Jk 

9 

o 


o 

o 

9 

O 

o 

© 

(O < 

r-| 

fn 

• 

• 

• 

• 

ft 

• 

ft 

• 

ft 

ft 

ft 


ft 

ft 

ft 

ft 

ft 

ft 

ft 

-ft 

• 

ft 

ft 


1 


««« 



««4 

«>« 

o 

o 

9 

9 

9 

9 

CM 


CM 

/M 

•ft 

O 

9 

o 

9 

9 

9 

9 


« 











a 

















X) 

>o 

X) 

CM 


o 

f- 

■n 

TM 

9 

9 


9 

-9 

9 

n 

•ft 

•ft 

•ft 

•ft 

•ft 

•ft 

9 


g 


• 

• 

♦ 

• 

ft 

ft 

• 

ft 

.ft 

• 

ft 

CM 

• 

• 

ft 

ft 

• 

ft 

ft 

ft 

ft 

-ft 

4 




wof 





ftft 

o 

9 

9 

9 

9 

-M 

"A 

CM 

•ft 

o 

9 

•9 

o 

9 

© 

o 

9 


in 



















a 

1 

• 

1 

1 



'• 




o 

kn 

. 

•Ift 

9k 

n 

*n 

'9 

O 


in 

o 

M- 

9 

9 

CM 

•ft 

•ft 

•ft 

•ft 

O 


m 


• 

• 

« 

ft 

^ ft 

'ft 

-ft 


ft 

ft 

-ft 

—4 

ft 

ft 

ft 



ft 

ft 



ft 

ft 


II 


fM 


fM 


•>t 


o 

9 

© 

O 

9 


•ft 

••4 

9 

O 

9 

9 

o 

O 

9 

9 

o 




















f 

1 

a 

t 

a 



W 

X 



9 

««4 

9» 

■» 

- 

o 

x> 


Ml 

*M 

•ft 


- 

Cl 

Ok 

9 

- 

9 

9 

M* 

*A 

CM 

- 


131 



Table 2A-12(c) Vertical Wind Speed. W^ 

Case 12 (L) Series No. 1759, 

Length of Field L = 3636.0 nii, Horizontal Increment Ax = 90.90 meters 
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ORIGINAL PAGE IS 
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Table 2A-13(b) Wind Speed Perpendicular to Frontal Motion,Wy 

Case 13 (M) Series No. 0211, 

s”^. Length of Field L = 6349.2 m. Horizontal Increment Ax * 158.73' meters 
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Table 2A-13 (g) Vertical Wind Speed, 

Case 13 (M) Series No- 0211, 

Length of Field L = 6349.2 m. Horizontal Increment Ax = 158.73 meters 


of page is 
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Table 2A-14(a) Wind Speed in Direction of Frontal Motion, W^ 

Case 14(N) Series No. 0109, 

8.0 m s~^ Length of Field L = 5128.4 m , Horizontal Increment Ax *= 128.21 meters 
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Table 2A-14(b) Wind Speed Perpendicular to Frontal Motion, Wy 
^ Case 14 |N) Scries No. 0109, 

8.0 m s , Length of Field L = 5128.4 m., Horizontal Increment Ax = 128.21 meters 


ORIGINAL PAG^ IS 
OF POOR QUAUTX 
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Table 2A-T5(a) Wind Speed in Direction of Frontal Motion, 

Case 15(0) Series No. 0436, 

11.4 m s \ Length of Field L = 7407.6 m. Horizontal Increment Ax = 185.19 meters 
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Table 2A-15(b) Wind Speed Perpendicular to Frontal Motion, Wy 

Case 15(0) Series No. 0436, 

11.4 m s ^ Length of Field L = 7407.6 m , Horizontal Increment Ax = 185.19 meters 
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Table 2A-15(c) VertiGal Wind Speed, 

Case 15(0) Series No. 0436, 

11.4 m s , Length of Field L = 7407.6 m , Horizontal Increment Ax = 185.19 meters 
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Table 2A-16(a) Wind Speed in Direction of Frontal Motion, 

Case 16 (P) Series No. 0220, 

6.9m s \ Length of Field L = 4444.4 m. Horizontal Increment Ax = 111.11 mete: 


•* 


01 

U 
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Table 2A-16(b) Wind Speed Perpendicular to Frontal Motion, 

Case 16 (P) Series No. 0220, 

6.9 m s~^. Length of Field L = 4444.4 m. Horizontal Increment Ax = 111.11 meters 
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Table 2A-16(c) Vertical Wind Speed, W^ 

Case 16 (P) Series No. 0220, 

6.9 m s~^. Length of Field: L = 4444.4 m , Horizontal Increment Ax = 111.11 meters 
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Table 2A-17(a) Wind Speed in Direction of Frontal Motion, 

Case 17 (Q) Series No. 1837, 

ns , Length of Field L = 11428.4 m , Horizontal Increment Ax = 285.71 meters 
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Table 2A-17(b) Wind Speed Perpendicular to Frontal Motion, W^ 

Case 17 (Q) Series No. 1837, 

= 17.5 m s ^ Length of Field L = 11428.4 mv. Horizontal Increment Ax = 285.71 meters 
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Table 2A-17(c) Vertical Wind Speed, 

Case 17(0) Series No. 1837, 

-1 

17.5 m s , Length of Field L * 11428.4-^, Horizontal Increment Ax = 285.71 meters 
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Table 2A-18(a) Wind Speed in Direction of Frontal Motion, 

^ Case 18 (R) Series No. 1651, 

25.0m s , Length of Field L = 16000.0 m. Horizontal Increment Ax = 400.0 meters 



148 



Table 2A-T8(b) Wind Speed Perpendicular to Frontal Motion, W^ 

Case 18 (R) Series No. 1651, 

-1 . 

25.0 m s , Length of Field L => 16000.0 m. Horizontal Increment Ax = 400.0 meters 
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Table 2A-19(a) Wind Speed Jn Direction of Frontal Motion, W^ 

Case 19(S) Series No. 1904, 

11.0 m s Length of Field L = 7142.8 m. Horizontal Increment Ax = 178.57 meters 
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Table 2A-19(b) Wind Speed Perpendicular to Frontal Motion, Wy 

Case 19 (S) Series No. 1904, 

11.0 m s Length of Field L = 7142.8 m , Horizontal Increment Ax = 178.57 meters 
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Table 2A-19(c) Vertical Wind Speed, 

Case 19 (S) Series No. 1904, 

11.0 m s~^. Length of Field L = 7142.8 m, Horizontal Increment Ax = 178.57 meters 
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Table 2A-20(a) Wind Speed in Direction of Frontal Motion, 

^ Case 20 (,T) Series No. 0029. 

19.9 ms. Length of Field L = 13333.2 m. Horizontal Increment Ax = 333.33 meters 
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Table 2A-20(b) Wind Speed Perpendicular to Frontal Motion 
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APPENDIX 2B 

GRAPHICAL ILLUSTRATION OF THUNDERSTORM WIND SPEED PROFILES 

This appendix contains graphical illustrations of the wind speed pro- 
files encountered along flight paths having a 3° glide slope through each of 
the 20 thunderstorm cases. The position of each flight path relative to the 
streamlines fixed in the frame of reference moving with the storm is illus- 
trated in each of the figures, and is chosen to terminate in the left-hand 
corner of each data set. The flight paths are, therefore, essentially arbi- 
trary flight paths which could be encountered during any routine landing 
through a passing thunderstorm. 

The longitudinal, lateral, and vertical wind speeds encountered along 
the glide slope are shown. The ordinate in the wind speed profiles is 
height, z, nondimensionalized with the length, H = L tan 3°, where L is the 
length of the wind field. Each profile is the wind as seen by an airplane 
traveling along the flight path drawn across the streamlines as shown in the 
upper figures. Note that the streamlines plotted are relative to the speed 
of the front which, for reference purposes, has been indicated on the hori- 
zontal wind speed profile with a vertical dashed line. The wind speed 
profiles are relative to the fixed earth frame of reference. 

The purpose of these illustrative profiles is to provide an impression 
of the diversity and severity of wind shear that can be encountered during a 
routine approach when a thunderstorm is in the vicinity of the terminal area. 
The data sets do not necessarily represent a wind field through the center of 
the storm nor the worst downdraft conditions, but represent storms which may 
have passed either to one side or the other of the tower from which the data 
were measured. Therefore, the profiles do not illustrate necessarily the 20 
most severe flight paths which would be encountered through the given thun- 
derstorms but rather a collection of wind speed profiles that represent an 
averaged situation. 
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Figure 2B-3 Typical Wind Profiles along a T Glide Slope: 
W = 8.6 m s"\ H 5714 tan 3“ ("0 
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Figure 2B-5 ^pical Wind Profiles along a 3* Glide Slope: Case 5(E) Series No. 1924, 31 May 71, 

VJ = 16.7 m s"^, H = 11428 tan 3* (m) 
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Figure 2B-8 Typical Wind Profiles along a 3° Glide Slope: Case 8(H) Series No. 1712, 23 May 74, 

W = 8.5 m s"\ H = 5556 tan 3* (’>") 



















2B-1 










169 








170 


Figur 







171 


2B-1 

















173 




tOUd pu 





















176 








ORIGINAL PAGE IS 
OF POOR QUALITY 




appendix 2C 

COMPUTER PROGRAM FOR CALCULATION OF TURBULENT 
THUNDERSTORM WIND FIELDS 

This appendix describes a computer program used to obtain the longi- 
tudinal, lateral, and vertical wind speeds at a given point (x,z) and all six 
velocity gradients at that point. The computer program has the option of 
superimposing a non-Gaussian turbulent field on the thunderstorm wind speeds. 
When the turbulence option is utilized, an input for time is required. The 
specified time designates to the program the period over which the turbulent 
fluctuations are computed. 

The subroutine VELO calls seven additional subroutines as schematically 
illustrated in Figure 2C-1. These seven subroutines are entitled TURB, 
ROMDON, FILTER, GAUSS, FFT, WIND, and BUILD. 

Program Description 

The program is written in FORTRAN IV and generates the wind field with 
turbulence superimposed for (x,z,T) inputs from the calling flight simulator 
program. 

Subroutine VELO 

For the given position and time the VELO subroutine furnish the wind 
velocity [U,V,W output ] and velocity gradient [UDX,UDZ,VDX,VDZ,WDX,WDZ, 
output ] in a flow field which is one of 20 cases of cold air outflow from 
thunderstorms detected by National Severe Storms Laboratory during 1971 
through 1974 [2-1]. 

Also, a non-Gaussian Dryden spectrum turbulence fluctuation is added to 
the wind velocity. 

VELO first calls the subroutine FILTER which sets up the transfer func- 
tion for tailoring the random signal. It then calls subroutine TURB. The 
subroutine TURB processes the random signal generated in subroutine ROMDON 
through the prescribed filter system and outputs fluctuations w^, w^ and w^ 
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Figure 2C-1 Schematic of VELO Subroutine 





having a non-Gaussian distribution with a Dryden spectrum form (Equation 2-4). 
The transfer to VELO is achieved with a common statement. VELO then calls 
subroutine WIND to obtain the mean wind velocities at (x,z) to which the 
corresponding fluctuations are added. The wind speed and wind speed 
gradients are then computed and the data is output in units of m s"^ and s"\, 
respectively. 


Nomenclature 

Subroutine VELO(X.Z,T,HAM,NBB,U,V,W,UDX,UDZ,VDX,VOZ,WDX,WDZ) 

X Position vector in longitudinal direction (user's units) [INPUT] 

Z Position vector in vertical direction (user's units) [INPUT] 

T Controls time period of turbulent signal (sec) [INPUT] 

HAM Ratio of unit (meter/user's units) [INPUT] 

NBB I Integer number for calling a new thunderstorm case. [INPUT] 

NBB=1 calls a new thunderstorm data set into storage [OUTPUT] 

U Velocity in longitudinal direction (m s"^) [OUTPUT] 

V Velocity in vertical direction (m s“^) [OUTPUT] 

W Velocity in lateral direction (m s”^) [OUTPUT] 

UDX Longitudinal velocity gradient in x direction (s”^) [OUTPUT] 

UDZ Longitudinal velocity gradient in z direction (s'^) [OUTPUT] 

VDX Vertical velocity gradient in x direction (s"^ ) [OUTPUT] 

VDZ Vertical velocity gradient in z direction (s"^) [OUTPUT] 

WDX Lateral velocity gradient in x direction (s~^) [OUTPUT] 

WDZ Lateral velocity gradient in z direction (s"^) [OUTPUT] 


T is incremented with each integration time step of the calling program 
and is a control variable for the length of the random, turbulent signal 
generated. When VELO is first called in it generates a turbulent signal at 
height z for a specified time period TEND. Each time VELO is called follow- 
ing the initial case, i.e., NBB = 1, T is compared with TEND, when T ^ TEND a 
new turbulent signal at the current height z of time period TEND is generated. 
Thus a quasi-turbulence variation with height is achieved. 

NBB is set equal to unity for the first time VELO is called. It auto- 
matically changes NBB = 2 until another thunderstorm is desired at which time 

the calling program must reset NBB = 1. 
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Listing of Subroutine VELO 


ORIGINAL PAGE IS 
OP POOR QUALErai 


aUBRuUTiNb Vt:LOCX,Z, IV, HAM,I'JBB,U,V,W,UUX,UDZ, VDX, VOi:,, wDX, WOZ) 

CGMMON/Sr2/IX 

C0rtWUlM/SA'4/HC3,3, 12«) 

CO«MOiM/ST7/DX 
COMMyN/SX12/DX,NNN 
COMMON/ST13/GGC3, 128) 

IK(NBB-l) 101,101,105 
101 lX=b5549 

CAuU KIuTBRCZ) 

CALL rUKB(Z) 
rti4U=0. 

NCC = 0 

103 i't;NO='f Ei'ID + OT^RNn 

105 XFCr-TtRl)) 109,108,108 . 

108 call TUHB(Z) 

RCC=LCC+i 

GO TO 103 

109 tH=T/DT 
.Mr=eR+i-NiMi\i*i\icc 
U t>=G G (l , N X ) 

VP=GG(2,i'ri') 

rtP=GGC3,ivT) 

CALL WIND (X , Z , HAM , MSB , WX , WZ , w Y , WXX , wXZ , wZX , wZZ , W YX , R Y Z ) 

J=WX+UP 

V = WZ + \?P 

'-H = WY + WP 

NU=i'JT + l 

UDX=WXX+(GG( 1 ,nU)-UP)/DX 

VDX=RZX+(GG(2,NU)-VP)/DX 

vKDX = WYX + (GG( 2 ,tMU)-«fllP)/DX 

UDZ=WXZ 

YUZ=WZZ 

wOZ=WYZ 

REXUKN 

END 


Subroutine TURB 

The subroutine TURB is called by the subroutine VELO. This sub- 
routine carries out the calculation for the turbulent fluctuations based 
on non-Gaussian turbulence model of Reeves, et al . [2-3]. A Dryden 
spectrum is used. A schematic of the non-Gaussian turbulence simulation 
model is given in Figure 2C-2. 
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Figure 2C-2 Physical Interpretation of a Single Component of 
the Non-Gaussian Turbulence Model 


TURB first calls the subroutine ROMDON to obtain three independent dis- 
crete Gaussian white noise random functions N^> N|^, N^. These signals are 
processed through a common statement by the subroutine FILTER having the 
filter function 

After mathematical analysis the simulated turbulent fluctuation w (t), 

A 

w (t) and w (t) are determined and with a common statement returned to VELO. 

y z 

Listing of Subroutine TURB 
dUdKOUrirtlfc TURd(Z) 

DlMfiNSlON RRCi, 128) ,R1(3/128 J ,XRU2B) ,XU i2b) , c.,H ( 3 , 3 , 1 2 b ) 

COMMON/SX2/1X 

COMMON/Sf4/H(i,3, 128) 

CU.4MUN/StH/NT2 
CDMMON/S'n2/D'r,NNN 
CUrt(^ON/ST13/GG(3, 128) 

«r=NNN 

R=l, 

IT=1 
99 J = 1 

CALL RUMOUM(RR,RI,Nr,Nr2) 

100 00 110 K = l,N'f 

XR(K)=H( IT, J,K)*RR( J,K) 

XI CK ^ = Q. 

110 cuNriNUfc 

call FFT(XH,X1,WT,NT2/2) 

DO 111 K=1,NT 
GriUT, J,K)=XH(K) 

111 continue 
J=J + 1 

lE(J-3) 100,100,130 
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ORIGINAL PAGE IS 

130 DO 140 KS1,NT OF POOR QUALITY 

XR(KJ=GHUT,1 rK)*GH(II,2,K) 

GG(IT,K)=K/( l.+rt»R)»*0,5*XR(K)-*-l ./(l.+H^Rl^^O.b^GHdTrS/K ) 

140 COrtTlNUt 
XI=1T+1 

IF(lT-3) 99,99,150 
150 CONTINUE 
«RITE(b,2) 

DU IbO 1=1,3 

WRITECb, 1) (GG( 1 ,K) , K = 1 , t\ll ) 
lF(i-2)lb4, 155,lb0 

154 rtKIT£(b,3) 

GO TO IbO 

155 WRITE(b,4) 

160 CONTINUE 

1 FORMAT CIO (IX, EU. 4)) 

2 F0RMATC15X, 'LONGITUlNAt,' ) 

3 FURMAT(15X,'VERTlCAt.») 

4 FORMATClbX, 'LATERAL' ) 

RETURN 

END 


Subroutine ROMDON 


The subroutine ROMDON is called by the subroutine TURB. This subroutine 

calls the subroutine GAUSS which generates a random Gaussian white noise 

signal. Three series of normalized zero mean random function Ng(t), N^^Ct) 

and N (t) are transferred to the subroutine FFT which fast Fourier transforms 
c 

N (t), N. (t) and N (t) and returns the transformed signal to TURB. 

Si D C 


Nomenclature 

Subroutine ROMDON (RR,RI.NT,NT2) 

RR(1,J) Real part of Fourier transform of Ng(t) [OUTPUt] 

RR(2,J) Real part of Fourier transform of N|^(t) [OUTPUT] 

RR{3,J) Real part of Fourier transform of N^(t) [OUTPUT] 

RI(1,J) Imaginary part of Fourier transform of Ng^(t) [OUTPUT] 

RI(2,J) Imaginary part of Fourier transform of Nj^(t) [OUTPUT] 

RI(3,J) Imaginary part of Fourier transform of N^(t) [OUTPUT] 

NT Integer number of discrete noise signal, [INPUT] 

usually NT = 128 


NT2 (2)*'*^^ = NT (for example, if NT = 128, NT2 = 7) 
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Listing of Subroutine ROMDON 


SU^KOUXlNt, RUMOnN(RR,Rl,NT,NT2) 

COMMON/SX2/XX 

UXrtfc..<iSION KR(3,12B),RI(3,1^B) , XR ( 1 2B ) , XX Cl 2B ) 

00 110 1=1,3 

SUM=0. 

00 105 0=1, NX 

CALL GAUSL(IX,1. ,0.,R) 

RHCI,JJ=R 
SUM=SUM+HK( X ,0 ) 

10 5 C UiM l'XN Ut 
AVt;=SUM/NX 

sxu=o. 

DO 106 0=1, NX 
RK(l,JJ = KK(l,0)-AVii: 

SIG=SXG+RR(I ,0)*»2 

106 CONXXNUt 
SXG=(SIG/NT)*=»=0.5 
00 107 0=1, NX 
HR(X, J)=Hk(X,0)/SXG 

107 CONXlNOc: 

OQ 108 0=1, NX 
XK(0)=KH(X, 0) 

X1CJ)=0. 

108 CONXINUC 

CALL KFT(XR,XI,NT,NX2,n 
DO 109 0=1 , NT 
KR(X,J)=XK(0) 

HI(l,J)=Xi(0) 

109 CONXlNOe 

110 CONXINUfc: 

RETURN 

END 


Subroutine GAUSS 

The subroutine GAUSS is called by the subroutine ROMDON. Each time 
GAUSS is called it generates a new random number. Random numbers generated 
by GAUSS are white noise having a Gaussian distribution. 

Nomenclature 


Subroutine GAUSS(IX.S,AM,V) 

IX Arbitrary starting number, IX = 65549 

S Adjusting number, setting S = 1.0 gives standard 
- deviation of unity 


[INPUT' 

[output; 

[INPUT] 
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AM Adjusting number, setting AM = 0.0 gives zero mean 
V Random number output 


[INPUT] 

[OUTPUT] 


Listing of Subroutine GAUSS 


SUBRUUTI Wfc GAUSS ( IX , S , AM , V ) 

REALMS X 

A=0.0 

DO 50 1=1,12 
1I=IX*65539 
XFUY) 5,b,6 

5 iY=iy+2i474b3b47+l 

6 Y=IY 

Y=Y*0.4bb66l3D-y 
1X^1 Y 
50 A=A+Y 

V=CA-b,U)*6+AM 

KErUHN 

ENU 




Subroutine Filter 

The subroutine FILTER is called by the subroutine VELO. The input list 

contains the height z. This subroutine sets up nine filter functions shown 

in Table 2C-1, for generating the three direction fluctuating signal, w , w 

X y 

and w^. 

These are returned to TURB through a common statement where H(IT,0,K) 
equals H^CK) if J = 1, H. (K) if J = 2 and H_(K) if J = 3. The value of 

a P C 

IT = 1 specifies longitudinal fluctuations; IT = 2 vertical and IT = 3 
lateral. 

Listing of Subroutine FILTER 

SUBROUri'\iE FILTER (Z) 

COMMON/St4/H(3,3,128) 

COMMON/ST12/DT,NjyiN 

NT=NNN 

Zl=Z*3,2ai 

UMEAN=8. 

ALUs ( Zl/ ( 0. 177+0. 832*0. 001 ♦Zl)**l. 2) /3. 281 
A L Vs ALU 
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TABLE 2C-1 


TRANSFER FUNCTIONS OF THE NON-GAUSSIAN MODEL 
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ALrtsALU 

AMUs 0 . 1 0 &♦ 8 . / C 0 . 1 7 7 + 0 . S3 2 ♦ 0 . 0 01*Zl ) **0 . 4 


amv=amu 

AMWsAMU 

UOUsALU/UMtAN 

VOUsaLV/UMCAN 


WOUsAliW/UMEAN 

ORIGINAL PAGE IS 

0D=12«»*0.5 

OF POOR QUALITY 

0£s3.**0.t> 
FC=1 ,/(2.*UT) 
DO 110 K=1,NT 
SsKC*(K-U/MT 


rtCl , 1 ,K)=4.*AMU*U0U/(1 .+2.»UUU*S) 
H(l,2,K)=l./(l.+2.*UOU*S) 
H(l,3,K)=AMU»(2.*UOU)=i‘*0.5/(l. + UtjU»S) 

H(2, 1 ,K)=AMU*D0*V0U**2/( 1 .+2.*V0U»S) 
H(2,2,K)=S/(l. + 2.»VQJ*S)»=i'2 

H(2, 3,K) = AMV*VOU**0,5*(1 ,+13fc:*V0U*£>J/(l . + tfUU*S)*^2 
HO, l^K)=AMW*DD»WOU**2/(l. + 2.*WOU*S) 
h(3,2»K)=S/(1.+2.*WOU*S)#*2 

H(3,3rK)=AMw*WOU*«0.5»(l.+Dti;*W(JU»5)/(l, + wUU^S)*»2 
110 CONTINUfc; 

RCXURN 

£NU 


Subroutine FFT 

The siibroutine FFT provides a discrete fast Fourier or inverse Fourier 
transformation and it is called by the subroutine ROMDON and the subroutine 
TURB. 

In the subroutine ROMDON, FFT is utilized to transform and 

from the time domain to the frequency domain. 

In subroutine TURB, after the frequency function has passed through the 
filter function, FFT is used to transform the frequency function back to the 
time functions a(t), b(t) and c(t). 

Nomenci ature 

Subroutine FFT(XR,XI,N,NU,NGG) 

XR Array of the real part of the function to be transformed 

XI Array of the imaginary part of the function to be 

transformed 

N Total number of the Array (N); usually N=128 
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[INPUT' 

[OUTPUT' 

[INPUT' 

[output; 

[INPUT] 



NU (2)NU=N; usually NU=7 which gives N-128 

NGG If NGG=1, the subroutine performs the Fourier Transform 

If NGG=2, the subroutine performs the Inverse Fourier 
Transform 


Listing of Subroutine FFT 


SUBROUTllfK FFT(XR,XI,N,NU,NGG) 
UXMtNSION XK(N),XI(N) 
IFCNGG.E^J.n GO TO 321 
00 32 0 1 = 1, iM 

320 XI(n=-Xl(I) 

321 N2=N/2 
NU1=NU-1 
K=0 

00 lod L=1,NU 

102 DO 101 1=1,IM2 
P=lblT8CK/2»*NUl ,N0) 

AHG=6.28 31H*j*P/f UOAKN) 
C=COS(ARG) 

S=S1N(ARG) 

M = K+1 
K1R2 = M + R2 

rH=XR(KlN2)#C+XI(KlN2)*S 
Ii=XI(KlN2)*C-XR(KlN2)*S 
XR(K1N2)=XRCK1)-TR 
XUK1W2)=XI(K1)-TI 
XKCKl )=XR(K1 )+rH 
XI(KU=XXCK1)+T1 
101 K=K+1 
K=K+N2 

IFCK.UX.iO GO TO 102 
K = 0 

NJ1=NU1-1 
100 N2 = N2/'2 

DO 103 K=1 ,N 
I=I8ITR(K-1 ,NU)+1 
IF(l.LE.K) GO TO 103 
TR=XRCK) 

II=X1(M 

XRCK)=XR(1) 

Xl(KJ=XItl) 

XRC1)=TR 

X1(X)=IX 

103 CONIMnUE 
IF(NGG.EU.I) go to 121 
DO 124 X=1,N 
XRCX)=XR(X)/N 
XXCX)=XXC1)/N 


[INPUT] 

[INPUT] 
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J24 CONTINUE 
121 RETUKN 
ENO 

FUNCTION 1B1Tk(J,NU) 

JlsJ 

IBITHsO 

DO 200 1=1,NU 

J2=dl/2 

IbllR=l6J.lK»2 + (Jl-2*J2} 
200 01=02 
KEXURW 
END 



Subroutine WIND 

The WIND subroutine is cal Ted by the subroutine VELO for the given 
position (x,z) from VELO. Interpolation by an area-weighting method for the 
velocity at the point x,z in the flow field based on the tabulated data for 
the given thunderstorm data set placed in storage by subroutine BUILD. The 
area-weighting interpolation method is illustrated in Figure 2C-3. 

The area-weighting scheme calculates the velocity and velocity gradient 
of each point by using the four nearest neighboring grid point velocities. 



Figure 2C-3 Area-Weighting Technique 
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The velocity at point P is given by: 
where 


A = A-j + A2 + Ag + 

The same interpolation method is used for the velocity gradients. 

Nomenclature 

Subroutine WIND (X,Z,HAM,KCK,WX,WZ.WY,WXX,WXZ,WZX,WZZ,WYX.WYZ) 

(X.Z) Position vector (users units) [INPUT] 

HAM Ratio of units (users units) [INPUT] 

KCK Integer number for storing a new thunderstorm case. [INPUT] 

If KCK=1 , the subroutine BUILD is called and a new 
thunderstorm data set is read in; otherwise KCK=2. 

WX Velocity in longitudinal direction 

WZ Velocity in vertical direction [OUTPUT] 

WY Velocity in lateral direction 

WXX,WXZ Velocity gradient 

WZX.WZZ Velocity gradient [OUTPUT] 

WYX.WYZ Velocity gradient 


Listing of Subroutine WIND 


SUbAUUriNt; WlND(X,Z,HAM,KCKf 

OlMENSiON DXX(2,2) ,DXZ(2/2) ,OZXC2,2) ,UZZ(2,2) ,UiX(2,2) , 
&Uf2C2,2) 

CJMrtUN/i’f/ft(41, n ,3) ,OX,DZ 

il f JKMaIC// , ' X-VALiUc: IN bUbKUUTlNb. iNfCR lfi» OUlSIDe FLOW 
S'KtiGIUN X = » 

21 FOKMaI (// , ' *♦* Y-VALUS IN SUtiRuUTiNS INIILR IS OUi’SIDt FLOw % 
S'KSGIUN K = ',ei3,b,//) 

31 FJKMAT(2X, ' ♦♦♦PUINI X = ',tl3.6, ' Z = St;i3.b,' IS OUXSIUE') 
IFCKCK.Nt.U LU TO 100 

call build 

KCK=2 
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100 

2 fc>=Z»HArt 

XF(XP.t,f.U..OH.ZP,lir.0,) GO iU 210 

OJ 101 1 = 2,41 

Al=Cl-l)»'JX 

llr(KP.Lfc;.Xn GU 10 102 

101 CJwlii>.ue ORIGINAL PAGE IS 

.KKlTt.C 6 , ll) XP OF POOR QUALITY 

102 1 P= 1-1 

OU 10 J J= 2 ,ll 
Zl=C 0 -l)* 0 Z 

IfXZP.Lt.Zll GU TO 104 
lOi CJNlii'iU£ 

wRlTfc:C 0 , 2 i| ZP 
104 JP=J -1 

ALiFA=(XP-UP- l ) *DX ) /DX 
Ot^f A = (ZP-( JP-U * 03 ) /DZ 
K = 1 

10 b VAL=( 1 . 0 -ALPA)*( 1 . 0 - 6 erA)*A(lP, JP,K)+tit.TA*( 1 . 0 -ALFA)»A(IP, JP 
$ + 1 ,K)+AUFA#( 1 . 0 -BKrA)*A(IP+l , JP , M + ALFA^Bfc'f A*A (IP+ 1 , JP + 1 , K) 
iF(K- 2 ) 10 b, 107 , 10 B 
10 b rtX=VAii 
10/ J»Z= )lkL 
K=K + 1 
GJ lU lOb 
lOti rt<=.VALi 

OJ 20B 1=1,2 
I UU 20« J=l,2 

lFCI,P,c, 0 ,l.AWt).I.EQ.l) Go 10 201 
iFUP.c. 0 . 40 .AlMU.I.EQ. 2 ) GU TO 202 

DXXCl, J)= 0 .b/ 0 X*(A(lP 4 l, JP+J- 1 , 1 )-A( 1 P+X- 2 , jp+o-i,in 
DZXCl, J)= 0 . 6 /UX*(A(lP+i, JP+J- 1 , 2 J-AUP+I- 2 , JP+J- 1 , 2 ) ) 

OlXti, J J= 0 .b/DX^( ACIPf I, JP+J- 1 , 3 )-ACXP+l- 2 , JP+J- 1 ,i)) 

GO TO 203 

201 L>XXCl,J) = tA( 2 ,JP + J-l,U-AU,OP+J-l,lJ)/OX 
DZXCl, J)=(A( 2 ,JP+J- 1 , 2 )-A(l , JP 4 J- 1 , 2 J)/ 0 X 
l)^XCl,JJ = (A( 2 ,JP + J-l, 3 )-A(i,JP+J-l, 3 J J/UX 
GO TO 203 

202 OXX( 2 ,J)s:CA( 4 i, JP + J- 1 , 1 )-AC 40 ,OP+J«- 1 , 1 ) J/OX 
UZXC 2 , J J=(A( 41 , JP+J- 1 , 2 )-A( 40 ,oP+J-l , 2 ) )/DX 
01 /Xt 2 ,d J = t A (41 , JP + J- 1 , 3 )-At 40 , JP+J- 1 ,31 J/OX 

203 XFCOP.EQ.i.AND.J.EQ.l) GO 10 204 
iFCJP.EQ. lO.ANO. J.E 1 . 2 ) G 3 lU 205 

OXZCI,d)=0.b/OZ*(A(IP+I-l ,dP+J, 1)-A(1P+1-1, JP+J-2,1) ) 

OZZCI, J)= 0 . 5 / 0 Z#(A( IP+X- 1 , JP+J, 2 )-A( 1 P+ 1 - 1 , JP+J- 2 , 2 )) 

U^Z C 1 -, d)= 0 . b/uZ*C AC IP + 1-1 ,JP+J, 3 ) -AC lP+l-l»JP + J- 2 , 3 ) ) 

GJ TO 20B 

204 UXZU,l) = CAClP + I-l»2,l)-ACiP + l-l,l,U)/0Z 
i)ZZCl,l)=CAClP+I-l,2,2)-ACiP+i-lrl»2))/UZ 
OiZCX, l)=CACiP+I-l ,2, 3)-ACiP+X-l , i ,3) )/DZ 
GJ 10 20 d 
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20b 0 X 4 U , 2 ) = ( A ( 1 P + X - 1 , 1 1 , 1) - A ( 1 P + 1 - 1 , 1 0 , 1 ) ) / D Z 
0 Z Z ( 1 , 2 ) = ( A ( 1 P + 1 - 1 , 1 1 , 2 ) - A 1 1 P + X - 1 , 1 0 , 2 J ) / 0 Z 
0xZ(l,2)=(A(lP+I-l , 1 1 ,3)-AtlP+I-l ,10, ij )/DZ 
20d CjiU’INUc: 

«AX=( 1 . 0 -ALFA)*( 1 . 0 - 0 ErA)*DXXCl,n+bElA»( 1 . 0 -AbFA)*DXX(l 

♦ ,2)+A0FA*t l.O-6tTA)*OXX(2, U + ALf A*tiEXA*0XXC2,2) 

WXZ=(l.U-ALf A)*(l,0-berA)»0XZa, l)+dt.lA»( 1 . 0 -ALfA)*DXZU 

♦ ,2)+AUt A^(1.0-Bb:rA)*0XZ(2, D + ALFA^bif A*DXZC2,2) 

rt ZX= ( 1 , U-At,FA ) ♦ ( 1 . 0-d£r A ) ♦OZX ( 1,1) +0E I’A* ( 1 . 0-ALF A ) ♦OZX ( 1 

♦,2)+AbFA^(1.0-8ETA) ♦0ZX(2,l)+ALFA*bt:iA»0ZX(2,2) 
^ZZ=(1.0-AuFA)*( l,0-Bei’A)*UZZCl, l)+ 8 El’A*C 1 . 0 -AbFA)*UZZ(l 

♦ ,2)+ALFA»a.U-bh.rA)»DZZ(2, 1) + ALF A»BfcXA*OZZ(2,2) 

.vlfX=(1.0-AL,f A)»( 1.0-BEf A)♦l>l^XClr U+bt.XA*U,0-ALFA)*6vXU 

♦ , 2)+AbFA»ll .U-BErA)»D^X(2, l) + AbFA*Bt*.XA«0XXC2, 2) 
rtxZ=(i.O-AbFA)*( l,0-6ErA)*DYZ(l , l)+BEl'A»a.O-ALFA)»OXZ(l 

♦, 2)+AbFA»(l .0-BETA)»UXZ(2, 1)+Abt A»BETA*DYZ(2,2) 
bU XU 21b ' 

210 W^lfE(6,il) XP,ZP 
21b REIUKN 
Ei'JO 

Subroutine BUILD 

The subroutine WIND is called by the subroutine BUILD which reads indi- 
vidual sets of wind speed data into storage. The thunderstorm cases are read 
in sequential order based on the numbering system utilized in this report; 
that is. Case 1 corresponding to Case A, Serial No^ 0446 of Reference [2-1] 
is read in first and remains in storage until KCK in list statement of sub- 
routine WIND is assigned the value 1. The subroutine WIND then calls the 
next thunderstorm case in numerical order. 

In subroutine BUILD, the grid system for the wind field is shown in 
Figure 2C-4 where IT=41 and JT=11. 



Figure 2C-4 Grid System of Wind Field 
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At point (I,J) the nortfial velocity W^(I,J) and vertical velocity Wj,(I,d) are 
stored as 

A(I,J.l) = Wj^(I,J); A(I,J,2) = W^(I.J) 

A common statement is used to transfer the data A(I,J,K) to the subroutine 


WIND. 


Listing of Subroutine BUILD 

Su:5KUUrii'i*i bUILl) 

Jii-lc.i'lSlUi'l ib(41) 

CJ-'iMu.M/rT/A(41 , 11 , 3) , l)X ,UZ 
1 e‘ JHMAf (20 a, 'HOK lZOMrAti V£L.QC11* 
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r/,/1 


i e Jit, viAX Cl 3Kb, 2) 

3 t J4iiAi’(4<ri2,2K,21F5. 1) 

4 t'jKiviAl'C/ »/ , 2UX, ' Vfc.RtlCAb VLLiUCm' 
b FJKMAr(2Ho.4,e‘b. 1) 

o KJRflAXCbA; 'uX =• ,F8.2,3X, 'DZ =' ,i-H.2,3X, ‘MtAM VX =Sf5.1) 
7 e JHMA1C3X, 'J/£* ,2X, 21(1.3, 2XJJ 
b KJkmAT(/,/, 20X, 'liAriiRAL M LhOCl'li' ,/ , / ) 

Kc.AL)(b,b) UX , DZ , UMtAi'i 
wKli‘t,(b,o) UX,OZ,UMb:AM 
i)j yo 1 = 1,41 
yu id(i)=i 
X = 1 


i'nKltKCo,!) 

lUU .VRXib(o,7) (lb( 1) , 1 = 1 ,41 , 21 
Oil 105 J = 1,U 

RiAJ(5,21 CA(I,J,K),l=l,4l) 

Lp‘C^,^Jt,2) <iU TU 10b 
UJ 101 1=1,41 
A(l, J,M=AU ,J,K)+U4fe.Ai'J 

101 CJi'll'lNOb 
10b CJNClNUb 

UJ 102 J=l,ll 
JaJ=12-J 

102 /iKllbCb, 31 JKJ,CA(I, J<J,K) ,1=1 ,41 ,2) 
|^ = ^ + l 

lfCK-3) 103,104,10b 

103 rtKl!b(o,4) 

fU 100 

104 w-ai’E(b,b) 
oJ l*U 100 

10b Rbl'UKN 
c,^u 
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APPENDIX 3A 


TABULATED STABLE AND NEUTRAL BOUNDARY LAYER DATA 

This appendix contains the tabulated data interpolated from the 
extensive measurement of wind speeds under varying conditions of stability 
reported by Clarke and Hess [3-1]. 

The tabulated values are the longitudinal and lateral wind speeds 

ZS/V/S. /S/V 

W (z) = W /u* and W^Cz) = W /u* for neutral and stable boundary layers for 
u > 0 where the dimensionless height z = zf/u*. The data actuany stored on 
computer cards and used with the computer lookup routine given in Appendix 3C 
i s 

AW^ = W^(z = 0.15) - W^(z) 
and 

AWy = Wy(z) - Wy U = 0.15^^ 

/N ^ 

These data have been converted to W and W using the relationships 

A y 

Wjz) = AWJ9 = 0.001) - AW fz) + W^{z = 0.001) 

X X XX 

and 

Wy(z) = AWy(2) - 4Wy(z = 0.001) 

where 

W^(z = 0.001) = K'^Un(0.001 Ro + 1) + 0.01125 y] 

X 

To extend the tables to values of z < 0.001, wind speeds can be computed from 

W (z) = — {£n(Ro z + 1) + 4.5 zy/K>; y ^ 0, 0 < z < 0.001 
X K 
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and 


ORIGINAL PAGE IS 
OP POOR QUALITY 


, /S * 

Wy(z) = 0 

The data are presented in a right-hand coordinate system with positive 
in the direction from left to right and W , positive into the plane of the 
paper. 

Tabulated results are given for values of p ranging from 0 to 200 in 
increments of Ap = 10 listed across the top of the tables and for values of 
ranging from 0.001 to 0.15 listed in the vertical left-hand column. 

From these tables the user can construct wind speed profiles for given 
values of p, u* and f. Note that z = zf/u* and Ro = u*/zf. 
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Table 3A-l(a) Longitudinal Wind Speed, W , Ro =10^ 
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Table 3A-l(b) Lateral Wind Speed, W^, Ro =10 
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Table 3A-2(a) Longitudinal Wind Speed, Ro “10^ 
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Table 3A-2(b) Lateral Wind Speed, Viy, Ro *10 
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Table 3A-3(b) Lateral Wind Speed, W^, Ro =10' 
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b 

-M, 

0 

-10.7 

- 1 0 . b 

-10. 

b 

-5.9 

-9,7 

-9,3 

-8 ,6 

-7.3 

-5.6 

0, 

.040 

-9. 

0 

-8. 

4 

-B.O 

-8.1 

-7. 

9 

-7,7 

-7. H 

-7.0 

-7.0 

-b.7 

-4.0 

u, 

. 0 io 

-6, 

3 

-5. 

7 

-b.2 

-5.0 

-b. 

0 

-5,3 

- b . b 

-b, b 

-5.5 

-4,2 

-2.4 

0, 

.020 

-4. 

0 

-3. 

1 

- 2 . 7 

-3,1 

-3. 

1 

-J.O 

-3,b 

-3.8 

-3.8 

-2.8 

-1.4 

0, 

.010 

-1. 

b 

-1. 

3 

-1,0 

-1.4 

-1 . 

4 

-1,2 

-l.b 

-1.8 

-1.9 

-1.4 

-0.6 

0, 

, OOb 

-0. 

8 

*•0. 

6 

-0.5 

- U . 6 

-0. 

0 

-0,5 

-0.7 

-0,8 

-0.8 

-0,6 

-0,3 

u, 

.001 

0. 

0 

0. 

0 

0.0 

0.0 

0. 

U 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


201 



A 4 

Table 3A-4(a) Longitudinal Wind Speed, Ro =10 



0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

U.1>0 

^4.b 

29.2 

3 3*9 

39.0 

4 2.9 

46*2 

4 b* 8 

51.5 

54,2 

5b. 5 

58.6 

0*140 

24. tt 

29.2 

33.9 

39.0 

42.9 

46.2 

48.8 

51.5 

54.2 

56.5 

58.6 

O.liO 

24. b 

29.2 

33.9 

39.0 

42.9 

46.2 

48.8 

51.5 

54.2 

56.5 

58.6 

0. r^o 

24.B 

29.2 

33.9 

J9.V 

4 2 .9 

46 * 2 

48.8 

51.5 

54.2 

56.5 

58*6 

o.uo 

24, b 

29 , 2 

33*9 

39.0 

42.9 

40.2 

4 b ,b 

51.5 

54.2 

56.5 

56*6 

0.100 

24 .b 

29*2 

33,9 

39.0 

42.9 

46 .2 

48.8 

51.5 

54.2 

56.4 

58.2 

0.090 

24. b 

29.2 

33.9 

39.0 

42,9 

4 o . 2 

48.6 

51,5 

54,2 

56.1 

57.6 

O.OtfU 

24.6 

‘29.2 

33,9 

39.0 

42.9 

4 b , 2 

48.0 

51.0 

53.0 

54.7 

56.2 

0.070 

24. b 

29.2 

33.9 

39.0 

42.6 

45.5 

47.7 

49.8 

51,7 

53,3 

S4.8 

O.OoO 

24.0 

29* 1 

33.8 

id. 8 

4 2.0 

44.1 

40.2 

48*2 

50.3 

52,0 

53*4 

0. oaO 

22.0 

26,7 

31.0 

35.9 

39.3 

42.1 

44.1 

46,2 

48,3 

50.0 

51.4 

0.040 

20.6 

24.5 

28.6 

33.3 

36.8 

39.6 

41.7 

43.7 

45*4 

46.8 

47.9 

O.OiO 

la.b 

22.5 

26.5 

30.9 

J4.1 

36.6 

3b*5 

40.2 

41.3 

42.1 

42.8 

0 . 020 

10.6 

20.5 

24.2 

27.7 

29.0 

29.2 

29 . 6 

29.8 

29,9 

31.0 

32.7 

0.010 

11 . b 

14,3 

16.5 

1 6.5 

19.0 

1 8.0 

18.6 

18.5 

18.4 

19.3 

20.6 

0. uob 

b.4 

9.7 

10.9 

11.9 

12*3 

12.3 

12.4 

12.5 

12*6 

13.2 

13.9 

0.001 

b.B 

6.0 

6.3 

0.0 

6.9 

7»2 

7.4 

7.7 

8.0 

8.3 

8.6 


100 

no 

1 2 0 

130 

140 

150 

160 

170 

1 80 

190 

- 200 - 

0.1 bO 

b8.0 

60 .0 

61,3 

62,2 

63. 1 

6 4.0 

64.3 

6 4.5 

64.8 

65*1 

65.4 

0.140 

68,6 

oO , 0 

61.3 

62.1 

62.9 

63.7 

63,8 

6 4.0 

64.3 

64,6 

64.8 

0.1 JO 

68.6 

oO. 0 

61.2 

ol .8 

6 2 , 6 

o3.2 

o3.3 

6 3.4 

6 3.7 

64,0 

64.3 

0. 120 

6 b, 6 

59.9 

60.9 

61.3 

61.9 

62.5 

62.6 

62.8 

63.1 

63.4 

63.7 

O.llO 

68.6 

59.6 

60.4 

60.7 

61.3 

62.0 

62.1 

6 2*2 

62.4 

62.6 

62.9 

0. loo 

68. 2 

69.2 

59.9 

6 0* 1 

o U . 7 

61.5 

61.6 

61,6 

61.7 

61.9 

62.2 

0.O90 

67.0 

58.5 

69.2 

69 . b 

6U.2 

60*9 

60.9 

60.9 

60.9 

61,1 

61.3 

O.OttO 

66.2 

67.3 

58.2 

5d. o 

59.2 

59.8 

69.9 

59 . 9 

59.9 

59 .9 

60.0 

0.070 

64. b 

5b* 0 

5b. 8 

57,2 

57,7 

58.4 

68.4 

58*4 

50.4 

58.5 

58.5 

0 . u o 0 

6 3.4 

54.5 

55.2 

65,4 

56.0 

50.9 

66.7 

So. 5 

56.6 

5o • b 

56.9 

0. ObU 

61.4 

62.6 

53.4 

53.5 

54.0 

64.8 

54.6 

54*4 

54.1 

54.0 

63.9 

0 . U40 

47,9 

4 9 .0 

49.6 

49 , O 

50.1 

bU. H 

50,0 

50.7 

50.5 

50 . 4 

50.4 

0.0 JO 

42.8 

43.7 

4 4*2 

4 3.9 

44*4 

46.1 

46 * 4 

4 5 . 0 

4 5.8 

4 o • 0 

46.3 

0,020 

J2.7 

34,0 

35*1 

35. 6 

36.4 

37.2 

37.5 

37.8 

38.0 

38.3 

38* 5 

0 . 01 0 

20.6 

21,7 

22.5 

23,2 

23.6 

24*6 

24*8 

25,0 

2 5*3 

2 5.6 

25.9 

0.00 b 

13.9 

14.6 

15.1 

15.5 

lo.O 

10.4 

le.7 

17.0 

17*3 

17.5 

17. 8 

0.001 

8.6 

8,8 

9.1 

9.4 

9,7 

1 0*0 

10.3 

10.5 

10.8 

11.1 

li.4 
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Table 3A-4(b) Lateral Wind Speed, Wy, Ro 





0 

10 

20 

30 

40 

50 

60 

70 

80 

90 

0 , 

. 150 

- 18.0 

- 20,4 

- 22.9 

- 25.6 

- 28.7 

- 32.0 

- 31.4 - 30.7 

- 29.8 

- 28.9 

0 , 

.140 

- 17.4 

- 19.8 

- 22.3 

- 24.9 

- 28.0 

- 31 .2 

- 30.6 

- 29.8 

-2 0.9 

- 27.9 

0 . 

.130 

- 17.1 

- 19.4 

- 21,7 

- 24.2 

- 27.2 

- 30.4 

- 29.7 ' 

- 28.9 

- 2«.0 

- 27.0 

0 , 

.120 

- 17.0 

- 19.1 

- 21.2 

- 23.4 

- 26.3 

- 29.5 

- 28.8 ■ 

- 28.1 

- 27.1 

- 26,1 

0 . 

.110 

- 16.3 

- 18.4 

- 20.5 

- 22.6 

- 25.3 

-28 . 3 

- 27.6 

- 26.7 

- 25.6 

- 24.6 

0 . 

.100 

- 15.5 

- 17.8 

- 19.9 

- 21.8 

- 24.3 

- 27.2 

- 26.2 - 25.2 

- 24.0 

- 23.0 

0 . 

.090 

- 14.5 

- 17.1 

- 19.4 

- 20.9 

-2j,2 

- 25.9 - 24.7 - 23.5 

- 22.3 

- 21.3 

u . 

Utf 0 

- 13.5 

- 15.5 

- 17.4 

- 18.9 

- 20.9 

- 23.2 - 22.3 - 21.4 « 

- 20.4 

- 19,6 

0 , 

,0 70 

- 12.1 

- 13.5 

- 15.0 

- 1 1>,5 

- lb . 6 

- 21.1 - 20,3 - 19.3 - 18.0 

- 17.2 

o< 

. OdO 

- 10.4 

- 11.3 

- 12.4 

- 13.9 

- 1 6 . 4 

- 19.4 - 18.4 - 17.2 - 15.2 

- 14.3 

0 , 

. 050 

- B .7 

- 9 .« 

- 10.9 

- 1 2 , u 

- 14,2 

- 17.0 - 15.9 - I 4 i 7 - 13.3 

- 12.3 

0 , 

, U 40 

- 7.0 

- 8,1 

- 9.2 

- 10.2 

- 12.1 

- 14.5 - 13.3 - 12.1 - 10.9 

- 9.9 

0 . 

,030 

- 5.3 

- 6.2 

- 7.2 

- 8.4 

- 10. 0 

- 11.6 - lu.b 

- 9.4 

- 8.2 

- 7.2 

0 . 

020 

- J . 7 

- 3.9 

- 4.4 

- 5.4 

- 7.1 

- 9.1 

- 7.9 

- 6.8 

- 5.8 

- 4.9 

0 . 

010 

- 1.8 

- 1.8 

- 2.0 

- 2.5 

- 3.6 

^ 4.8 

- 4,1 

- 3.4 

- 2.9 

“ 2.4 

0 . 

00 b 

- 0.8 

- 0.6 

- 0.9 

- 1.1 

- 1.6 

- 2.1 

-1 .« 

- 1.5 

- 1.3 

- 1.0 

0 . 

001 

0.0 

0.0 

0.0 

0.0 

0,0 

0.0 

OaO 

0.0 

0.0 

0.0 



100 

no 

120 

130 

140 

150 

16 U 

170 

180 

190 

0 , 

.150 

- 28. 0 

- 26.8 

-2 5.8 

- 25.2 

- 24.2 

- 23.0 

- 21 .8 

- 20.5 

- 19.1 

- 17.2 

0 , 

,140 

- 27.0 

- 25.8 

- 24.7 

- 24.2 

- 23.3 

- 22.2 

- 21.1 

- 19.8 

- 18.4 

-3 6.5 

0 , 

,130 

- 26.1 

- 24.8 

- 23.8 

- 23.3 

- 22.4 

- 21.4 

- 20.2 

- 19.0 

- 17.6 

- 15.7 

0 

.120 

- 25.2 

- 23.9 

- 22.9 

- 22.4 

- 21.6 

- 20.5 

- 19.3 

- 18.0 

- 16.6 

- 14.8 

0 . 

.110 

-2 3.7 

- 22.5 

-21 .5 

- 21.0 

-2 0 . U 

- 18.8 

- 17.9 

- 16.9 

- 15.6 

- 13.8 

0 , 

, 100 

-2 2.1 

- 21.0 

- 20.0 

- 19.6 

-1 0.0 

- 17. 3 

- 16.6 

- 15.8 

- 14.6 

- 13.0 

0 

.090 

- 20.4 

- 19.4 

- 18.5 

- lb . 1 

- 17.1 

- 15 - V 

- 15.3 

- 14.7 

- 13,8 

- 12.3 

0 

.080 

-l b . 8 

- 17.8 

- 17,1 

- 16.7 

- 15.8 

- 14.6 

- 14.0 

- 13.3 

- 12.5 

- 11.2 

0 

,0?0 

- 10.5 

- 15 . b 

- 15.3 

- 15,0 

- 14.2 

- 13.3 

- 12.7 

- 12.0 

- 11.3 

- 10.1 

0 

• UoO 

- 13.9 

- 13.5 

- 13.2 

- 13 ,0 

- 12.5 

- 11.9 

- 11.3 

- 10.7 

- 10.2 

- 8.9 

p 

.050 

- 11.5 

- 11.0 

- 10,7 

- 10.8 

- 10,5 

- 9.9 

- 9.7 

- 9.3 

- 8.6 

- 7.3 

0 

.040 

- 9.0 

- 8.4 

- 8.0 

- 8.1 

- V , 9 

- 7.7 

- 7 .b 

- 7,6 

- 7.0 

- 5,7 

0 

.030 

- 6.3 

- 5.7 

- 5.2 

- 5.0 

- 5.0 

- 5.3 

-5,0 

- 5.8 

- 5.5 

- 4.2 

0 

.020 

- 4.0 

- 3.1 

- 2.7 

- 3.1 

- 3.1 

- 3.0 

- 3.5 

- 3.8 

- 3.8 

- 2.8 

0 

.010 

- 1.8 

- 1.3 

- 1.0 

- 1.4 

- 1-4 

- 1.2 

- 1.6 

- l.b 

- 1.9 

- 1.4 

0 

.005 

" 0-8 

- 0.6 

- 0.5 

- 0.6 

- 0.6 

- 0,5 

- 0.7 

- 0.8 

- 0.8 

- 0.6 

0 

.001 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

O.U 

0.0 

O.O 

0,0 


100 

• 26 ,0 
- 27^0 
- 26.1 
- 25.2 
-? i .7 
- 22.1 
- 20.4 
- 10.8 
- 16.5 
- 13.9 
- 11.5 
- 9.0 
- 6.3 
- 4.0 
-l.« 
- 0.0 
0.0 

'200“' 

- 15.0 
- 14.3 
- 13 .6 
- 12.8 
- 11.0 
- 11.1 
-10. b 
- 9.0 
- 8.7 
- 7.2 
- 5.6 
- 4.0 
- 2.4 
- 1.4 
- 0.6 
- 0.3 
0.0 
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Table 3A-5(a) Longitudinal Wind Speed, Ro =10^ 



0 

10 

20 

30 

40 

bO 

fau 

70 

80 

90 

100 

0 . i :>o 

ly.o 

2 3.5 

28.2 

3 J .2 

37.1 

40.4 

43.1 

45.8 

48.4 

50.7 

52.8 

>^.140 

19.0 

23.5 

28.2 

33.2 

37.1 

40.4 

4 3.1 

45. tt 

48.4 

50.7 

52.8 

O.liU 

19.0 

?3.5 

2d. 2 

33.2 

37.1 

40.4 

43.1 

45.6 

48.4 

50.7 

52.8 

u.ii^O 

19.0 

2 3.5 

2e,2 

33.2 

37,1 ■ 

40.4 

43.1 

45.8 

48,4 

50.7 

52.8 

u.MO 

19.0 

23.5 

28.2 

33.2 

37.1 

40.4 

4 3.1 

45.8 

48.4 

50.7 

52.8 

U.lOO 

19.0 

23.5 

28.2 

33.2 

37.1 

40.4^ 

43.1 

45.8 

48.4 

50.6 

52.5 

u,uyo 

19.0 

23.5 

28.2 

33.2 

37.1 

yO.4 

43.1 

4 5.7 

4 8.4 

50.3 

51^9 

0, UdO 

19.0 

23.5 

28.2 

33.2 

37.1 

4 0.4 

42.9 

45.2 

47.3 

48.9 

50.4 

U.U7U 

19.0 

23.5 

28.2 

33.2 

3o,« 

39.7 

42.0 

4 4.1 

40. 0 

47.0 

49.1 

O.uou 

1«.9 

23.3 

28.0 

33.1 

36,2 

38.4 

40, 4 

42.4 

44.5 

40.2 

47.7 

0. J!>0 

lb. 9 

20.9 

25.2 

30.1 

33. b 

36.3 

3 b . 4 

40.4 

42.5 

4 4.2 

45.7 

0 . U *4 U 

14.9 

18.7 

22.9 

27.5 

31.0 

33.9 

36.0 

37.9 

39.7 

41.0 

42.1 

O.OiO 

12.« 

16.7 

20.8 

25.2 

2d. 3 

30.8 

32.8 

34.4 

35.5 

36.3 

37.0 

0.V20 

1 0 . b 

14.7 

1 8 .5 

2 1.9 

23.2 

23.4 

2 3. tf 

24.0 

24.1 

25.3 

26.9 

u.uiO 

b.O 

b.b 

10.8 

1 2 . 8 

13.2 

12.9 

12.8 

12.7 

12.6 

13.5 

14.9 

0. U05 

2.7 

3.9 

5.1 

0.1 

o.b 

0.5 

0.6 

0.7 

0.9 

7.4 

8.2 

O.uOl 

0.0 

0.3 

0.6 

0.8 

1.1 

1.4 

1.7 

2.0 

2.2 

2,5 

2.8 



It’O 

1 1 0 

120 

l3o 

140 

150 

160 

170 

1 tf 0 

190 

2 00 

0 . 15 U 

52.8 

54.3 

55.6 

5o,5 

57.3 

5b. 2 

5 b. 5 

58.8 

59.1 

59.3 

59.6 

U.14U 

52.8 

54.3 

55.5 

56.3 

57.1 

58,0 

5b.l 

58.2 

58,5 

58.8 

59.1 

0.130 

52, b 

54.2 

55,4 

56.1 

56.7 

57.4 

57.5 

57.7 

58.0 

58.2 

58.5 

0.120 

52.8 

54.1 

55.1 

55.0 

50. 1 

56.7 

5o, 9 

57,1 

57.3 

57,0 

57.9 

0.110 

52.8 

53.9 

54,7 

55,0 

55.5 

56.2 

50.4 

56.5 

56.7 

56.9 

57.2 

0.100 

52,5 

53.4 

54.1 

54,4 

55.0 

55.7 

55.8 

55,9 

55.9 

50.1 

56.4 

0.090 

51.9 

52.7 

53.5 

5 J. 8 

54.4 

55.2 

55,2 

55, 1 

55.1 

55.3 

55.6 

0.080 

50.4 

51.5 

52.4 

52,8 

53.4 

54.0 

54 . 1 

54.1 

54.1 

54.2 

54.2 

0.0 ;o 

49.1 

50,2 

51.1 

51.4 

52.0 

52, 7 

52.7 

52.7 

52.7 

52.7 

52.8 

0.0 oO 

47.7 

4 8.7 

49.5 

4:^ , 0 

50.2 

5i,l 

50.9 

50.8 

50.9 

51.0 

51.1 

U.U50 

45.7 

46.8 

47.6 

47.7 

48. 3 

49.1 

4 b . 8 

4 b. 6 

4 b. 4 

48.2 

48.1 

0.040 

42,1 

43.2 

43.9 

4 3, b 

44,3 

45.1 

4 5.0 

44. > 

44. b 

44.7 

44.6 

U.030 

37.0 

.37.9 

38.4 

38 , 2 

38.0 

39,4 

39.6 

39. b 

40.1 

4 0.3 

40.5 

0.020 

2o.y 

28.3 

29. 3 

29, y 

30 . 0 

31 .5 

31.7 

32.0 

32.3 

32.5 

32.8 

0.010 

14.9 

15,9 

16.8 

17,4 

1 b , 1 

lb, 7 

19.0 

19,3 

19.6 

19.8 

20 . 1 

0.005 

8.2 

8.8 

9.3 

9.8 

10.2 

10,7 

10.9 

11.2 

11.5 

n .8 

12.1 

O.OOl 

2.8 

3,1 

3.4 

3.7 

3.9 

4.2 

* 4.5 

4.8 

5.1 

5.3 

- 5.0 
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Table 3A-5(b) Lateral Wind Speed, W^, Ro =10^ 



0 

10 

20 

30 

4U 

50 

bO 

70 

80 

90 

100 

0.150 

-10.0 

-20.4 

-22.9 

- 2 5 . b 

-28.7 

-32,0 

- 3 i . 

4 

-30.7 

-29,8 

-28,9 

-28,0 

0.140 

-17.4 

-19.8 

-22,3 

-24.9 

-28.0 

-31.2 

-30. 

6 

-29.8 

-28,9 

-27.9 

-27.0 

U . i iO 

-17.1 

-19,4 

-21,7 

-24.2 

-27,2 

-30.4 

-29. 

7 

-28.9 

-28,0 

-27,0 

-26,1 

0 . 1 Zi) 

-17,0 

-19.1 

-21.2 

-23.4 

-20.3 

-29.5 

-28. 

8 

-2 b.l 

-27.1 

-26,1 

-25,2 

0.110 

-lb. 3 

-18.4 
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APPENDIX 3B 


GRAPHICAL ILLUSTRATION OF STABLE AND NEUTRAL BOUNDARY 
LAYER WIND SPEED PROFILES 

To provide an Illustration of the velocity profiles which can be 
encounter^ along a flight path during approach through neutral and stable 
boundary layers, seven velocity profiles have been computed using the model 
developed in this report. Both longitudinal and lateral wind speeds are 
shown. The conditions for which the velocity profiles have been computed are 
shown in Figure 3B-1 along with wind speed profiles calculated from the 
standard log-linear equation for the stable boundary layer, see References 
[3-2 and 3-5]. This enables a comparison of the simpler theoretical models 
with the more elaborate simulation models developed in this report to be 
made. One thing to note in comparing the velocity profiles is that the log- 
linear, theoretical models give only the longitudinal component of wind 
speed, whereas, the table lookup model provides both the longitudinal and the 
lateral components. The lateral wind speed component is quite significant in 
many of the cases illustrated and represents an estimate of the directional 
shear encountered during approach or landing. This directional shear cannot 
be obtained from the simpler theoretical equatidns. 

Figures 3B-2 and 3B-5 through 3B-7 are computed for a u* value of 
0.1 m s'^. This value of friction velocity when used to compute dimensional 
height from the nondimensional form z gives a value of z - 173 m at z = 0.15 
which is the maximum vertical extension of the tabulated data. Consequently, 
the aforementioned figures do not extend to the 500 m level as do Figures 
, 3B-1 , 3B-3 and 3B-4. 

Additional visualization of the velocity profiles which can be 
encountered in stable and neutral boundary layers with turbulence super- 
imposed are given in the text. Figures 3-8 through 3-13, Section 3.4. 
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Figure 3B-2 Stable Boundary Layer 



Figure 3B-3 Stable Boundary Layer 
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Figure 3B-4 Stable Boundary Layer 



Figure 3B-5 Stable Boundary Layer 
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Figure 3B-6 Stable Boundary Layer 
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Figure 33-8 Stable Boundary Layer 


APPENDIX 3C 


COMPUTER PROGRAM 

The subroutine STB is a FORTRAN computer program for calculating the 

wind speeds W and W and the wind speed gradients 3W /az and aW /az. Inputs 
X y X y 

to the program are height, z, friction velocity, u*, Coriolis parameter, f, 
stability parameter, p, and surface roughness, z^. These parameters must be 
introduced in units of meters and meters per second, respectively. The 
returned velocities and gradients have units of meters per second and inverse 
seconds, respectively. Wind speeds and gradients at a new value of height, 
z, are obtained by simply assigning a new value to z in the calling program. 

A new condition of stability is also achieved by assigning p a new value in 
the calling program. 

The user has the option of adding turbulent fluctuations to the wind 
speed by assigning the control integer NKK a value greater than two. 

Subroutine STB 

i 

The subroutine STB for a given height z and with the parameter, p, 

specified provides the wind velocity W and W . The user has the option of 

superimposing turbulence fluctuations on the mean wind speed if desired. The 

turbulence simulation routine uses the spectra developed in Section 3.5 and 

the z- transformation technique for generating the fluctuations w , w^ and w . 

X y z 

Subroutine STB first calls the subroutine INIT, which reads the input 
data into storage. The control variable NST is then set equal to two and 
further use of STB does not call the subroutine INIT. STB then calls the 
subroutine WINDF which interpolates wind speeds in the z direction for a 
prescribed value of p. 

At the user's option, NKK (0<NKK) can be set greater than two and the 
subroutine STB will call the subroutine DRYDEN which processes the random 
signal from the subroutine GAUSS and generates turbulent fluctuations having 
the modified Dryden spectra described in Section 3.4. These fluctuations 
are added to the turbulent wind field output. 


/ 
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Nomenclature 


Subrout 1 ne STB ( Z , AMU . USTR , F, ZO , NST ,WX . WY , DWX , DWY. DT , I X , NKK) 

1 Height (m) 

AMU Stability parameter (“333.33 < p < 216.67) 

USTR Surface friction velocity u* (m s"^) 

F Coriolis parameter, f (s"^) 

ZO Surface roughness, (m) 

NST Control variable NST = 1 initiates the storage of the 
tabulated wind speeds 

WX Wind speed in the x direction (m s“^) 

WY Wind speed in the y direction (m s”^) 

DWX Gradient of W^ in the z direction (s^^) 

DWY Gradient of in the z direction (s“^) 

DT Time increment At of the turbulent fluctuation (s) 

IX Initiating integer for random signal generator, value 

arbitrary 

NKK Integer for total terms of z transformation 
(if NKK ^ 2 turbulence is not added 
if NKK > 2 turbulence is superimposed on wind 
speed outputs) 


Listing of Subroutine STB 

' S U tt H U U T 1 .A) t S X ( Z P , A VI i J , u S X R , F , Z 0 , N I , W X , w ^ , U a X , U w X , D X 
C aMMUN/wXiVD/wOELC 34 > H , 2) , XMU , IMuM , AMUU ( 34 J ,UW(34 , 11 
$W<R(34) ,vvltR(34) , RAM J, VRfcF,UL)t,F 

COMMON /£>1/AL 
if( Ns r . Nt- . n cu TO 100 
IXsob54y 
CAUti INIT 
R0orK=U3i«/(F*Z0) 

N Si’s 2 

lOU COwTlNUti 
ZsZP+ZU 

CAUL sy I nuF C AMU , Z , WX , W X , OwX , DW Y , USTR , F , KUSTK , IZ) 
IFCnKK-2) 141,141,161 

141 OJsO, 

UV=0, 

CO TO 161 

161 CAUi, DHYUtN(Z,AMU,USrR,F, WX,DU,UV,Ort,DT,lX,NKK) 

Ibl wXsWX+UU 
W Y = W Y + OV 
HKTUKN 
CiMU 


[INPUT] 

[INPUT] 

[INPUT] 

[INPUT] 

tiNPUT] 

[INPUT] 

[OUTPUT] 

[OUTPUT] 

[OUTPUT] 

[OUTPUT] 

[INPUT] 

[INPUT] 

[INPUT] 


iX,NKX) 

2) f 
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Subroutine DRYDEN 


The subroutine DRYDEN first calls the subroutine GAUSS which generates a 
Gaussian white noise random signal. The random signal is then passed through 
a filter which generates the fluctuating output having a spectrum function of 
the Dryden form which fits the data of Kaimal [3-3]. The subroutine DRYDEN 
then returns the fluctuating wind components, DU, DV and DW to the subroutine 
STB. 

The subroutine DRYDEN utilizes the stability parameter, u, the friction 

velocity, u*, the Coriolis parameter, f, and the wind speed at height, z. 

2 

The scale frequencies, tIqj and turbulence intensities, a , are computed 
internally in the subroutine. 


Nomenclature 

Subroutine DRYDEN ( Z ,SMU , USTR , F , V , DU , DV , DW , DT , I X , NKK) 

Z Height (m) [INPUT] 

SMU Stability parameter, y [INPUT] 

USTR Friction velocity, u* (m s“^) [INPUT] 

F Coriolis parameter, f (s"^) [INPUT] 

DU Wind fluctuation in longitudinal direction (ms“^) [INPUT] 

DV Wind fluctuation in lateral direction (m s"^) [OUTPUT] 

DW Wind fluctuation in vertical direction (m s”^ ) [OUTPUT] 

DT Time increment of the turbulent fluctuations (s) [INPUT] 

IX Starting integer for generating random signal [SPACE] 

NKK Total number of terms for z transform [INPUT] 


Listing of Subroutine DRYDEN 


ai 

H2 
ai 
a 4 


SiJbK JU r L UK I Ofc.*v C/,,SMU,USTK,8,V,UU,DV,Dw,T,iX,NKK) 
xU0),V:CiUJ,K0(3j,i>lG(3),AKC3),AC(i) 
CDMMJi'i/Sr/ AI j 
ZAb = Z;/AL 

Kl=ZAb/( 1 ,+4.b*ZAL) 

'iF(ai-a. 028a)di ,fai ,a2 

Fuel) =0,01 4 

GO rj ai 


Fo( 1 )=o.t5>t=aj. 

iFca [-0.0176) a4,a4,ab 

F0(2)=0.0265 
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GvJ 1 0 H 6 

K0.( 2 ) =i .5*K1 

bo I r ( H I “ 0 , 0 3 4 i J c / , n / , » 6 
tj7 KO ( i ) = 0 , 0'9o2 
(j«.l f J H9 
8b f 0(3) =2. b ♦Hi 

89 iF(ZAb-i . ) ,9b ,yo,yo 

90 lc(2.Afj-1.22) 9/,y4,94 
9*4 .su;t j)=a. 

UfJ f'J 9 b 

9b SICiC 3) = 1 , 3-0. ib*ZAb^^0.b 
(in ivj 9 8 

9/ STC;(3)=6.49-b.32^ZAb 

9 b SI Q ( 2 ) = b i G ( 3 ) / (0 .1 / 7+0.002 7 / ♦ Z) ♦ ♦O . 4 
S.IG( I ) = biG(3) / (0.bb3 + 0.00139*ZJ4*D.b . 

O'.l 99 KK=l , 3 

A K ( K K ) = 8 b . 7 b ♦ V ♦ 1 0 t X K ) ♦ S X (G ( «< K ) ♦ 6 1 G ( K K ) / Z 
AC( KK) =2 3. 8b»\/»f 0 CxK)/Z 
99 Cib'jf lWUk; 

iOO 00 130 i = l,i'lXX 

Cnfjb GA JSS( iA , 1 . f 0. , H ) 

A( 1 )=H 

130 C^)3rii\aJF 

V ( 1 )= 0 . 

c=FKP(-Ac:(!Mbrt)*'n 

0= ( SDK r t A K ( .M db ) ) D ♦ c i . -C) / AC C Ortb ) 

141 1 = 2,i'!XK 

Y ( I J=C*K (i-U+0*XCl-l J 

141 Ctj.'] ri ivi'jit 

IcC-M3d-2) i42,l43»144 

142 = 

GO fO 143 

143 0 9=1 (iKK) 

Gi) I'O 14 b 

144 1)w=!i;(jKK) 

1 4 3 C I J O i' 1 l>i U t 

N n rt = (i rl S + I 

lK(MHb.bt.3) GU iu 100. 

Hfc,rUH.» 

CNO 


Subroutine GAUSS 

The subroutine GAUSS is called by the subroutine DRYDEN. Each time 
GAUSS is called it generates a new random number. Random numbers generated 
by GAUSS are white noise having a Gaussian distribution. 
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Nomenclature 


Subroutine GAUSS (IX,S,AM,V) 

IX Arbitrary starting number, IX = 65549 

S Adjusting number, setting S = 1.0 gives standard 
deviation of unity 

AM Adjusting number, setting AM = 0.0 gives zero mean 
V Random number output 

Listing of Subroutine GAUSS 


[INPUT] 

[OUTPUT] 

[INPUT] 

[INPUT] 

[OUTPUT] 


SOBKHUriMh: BAUS6(XX ,i>,A!'1, V ) 
.K£ALi*^i i 
A = 0 , 0 

DU 50 1=1, I Z 
I^ = U*bbg39 
1 K C i If ) 5 , b , b 
b i Y = i i +21 47 4b ibAZ+l 
h /slV 

k' = Y*') , Abbbbl 3 u-9 
£X = 1V 
b u A = A + y 

V=( A-b.O)^S+Aw 
HHJTiJRn 


origin A'Li 
OF POOR 


pA-GB 


Subroutine INIT 

The subroutine INIT is called once by the subroutine STB to initiate 
storage of the data. Subroutine INIT stores the data according to the grid 
system arrangement shown in Figure 3C-I. The data are stored as WDEL (I,J,K) 
where each row designated by J corresponds to the wind difference at that 
elevation as explained in Section 3.1. Each column designated by I repre- 
sents a prescribed value of y ranging from -333.34 to 216.67. The index K=1 
represents longitudinal wind speed and K=2 represents lateral wind speed. 

The stored data is transferred to the subroutine WINDF through a common 
statement. 
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2 3 1 34 


Figure 3C-1 Grid System: K=1 Corresponds to AW W (z = 0.15) - W^(z) 

and K=2 Corresponds to W = Q (z) - = 0.15) 

<j J J 



'J'j r I Me. iiMil 

(■ i V .1 I,v/ ,-J I .MU/ V7UC.L( i4 , I I . /.J , lMJ , IMUft , i 34 ) , Uw ( , 1 I , ■/ ) , 

Kk C 34 ) , M 34 ) * KAMU / ALil' A , UKbt , VKfcit , UUtf 
ZOO filR'AAHHti 10.4) 
aOO bnw.viA I'C 1 IF3. Z) 

i-c-c-U ifv. iiMlLAL DATA MU 
RFAi)(b,Z00) ( A.*iuu(l) ,i = l , 34) 
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IJU 2 K = l ,2 
DJ 2 1=1 , 34 

0 V ( i , 1 1 , K ) = C wubL 1 1 , 1 1 , K ) -1^01; L (1 , 10 , M ) / L)Z 
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RAMU=-400. 

AI.FA = 1 , 0 
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h 'j ‘ ) 
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Subroutind WINDF 


The WINDF subroutine is called by the subroutine STB for values of y and 
z. Interpolation is performed by an area-weighting method for the velocity 
at the elevation and for the stability condition, y, from the tabulated data 
placed in storage by subroutine INIT. The area-weighting interpolation 
method is illustrated in Figure 3C-2. 


okigwalpag^ 
OF POOR QUAUTY 


The area-weighting scheme calculates the velocity and .velocity gradient 
at the specified point by using the four nearest neighboring grid point 
values. The velocity at point P is given by: 

*^p " ^^1^1 ^ 2^2 ^3^3 

where 

A = A-| + A 2 + 

The same interpolation method is used for the velocity gradients. 

Nomenclature 

Subroutine WINDF(AMU,Z,WX,WY,DWX,DWY,USTR,F,ROSTR) 

AMU Atmospheric stability parameter, y [INPUT] 
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Figure 3C-2 Area-Weighting Technique 


z 

Height (m) 

Wind speed in the x direction (m s"^) 

[INPUT] 

wx 

[OUTPUT] 

WY 

DWX 

DWY 

Wind speed in the y direction (tn s"^) 
Wind speed gradient in z direction (s”^) 
Wind speed gradient in z direction (s"^) 

[OUTPUT] 

USTR 

Friction velocity, u*(m s"^) 
Coriolis paraineter (s"^) 

[INPUT] 

F 

[INPUT] 

ROSTR 

Rossby number, Ro 

[INPUT] 


Listing of Subroutine WINDF 

'■>JrW •,)(.• r 1 •'it AiHUir ( Al'iU, k. , A A , S.V Jf , UftX , IJW 1' , UlSl K , If ,KUSTH,lZ.) 

C .) »1 J I'i / A i .M l) / i'« 1 ; 1 1' C .5 4 , i l , / } f i ^ U I IhU W , M i»i U li C 3 A J , U IN t J A , 1 1 , id J , 
$ A.K ( 34) , W X lU 3‘0 , A>'iL) , Aor 4 > UKtt , VHtf’ , liiJb.F 
C ) "I .( • 'I / S r / A il 
i) M i,» = A i*i J •• fd A M U 
J.F( AiiSCDMU J -U , 01) 3 , A, 4 

4 1F( A'Ui. LI , ArtUM( 3 ) .(JK.A ■'lU.ijl, Ai»'ivJUC 34..) ) GO '1 0 W 
A L = 0 . 4 ^ J S T 4 / C r ^ A t‘i I.) ) 

K Ml 1 1 s A M 0 
1 = 1 

9 i=i+l 

1 1 *' ( A ' 1 0 « J r , AML) u ( I ) ) Go Ay y . 

1-1U=I 

I'^ U ■<'1 = 1-1 

A b KVm = A i) ’f 0 , 0 6 - 1 M 0 M t '/J 1 

li A r:F = AbF A^ nUKL C AMu ,1 rl J C 1 . - ALF A ) 4rtLitLL t iMUM » 1 » 1 ) 

VKFFsAtiF A^'NObb I J.MU , 1 , 3 J t ( 1-ALiif A) ’('WOtub 1 1 «UM , 1,2) 

•'(iK I rs C 6 , 1 6 ) UKFF , V KtiF 

lb fi IK '> A 1 ( 2 X , ' utd tc' = ' , 1 y . A , 3 A , ' V Htt*' =',Fy,4) 

1 1 .)K F = t A Li J G C KO 6 1 p * u .,-001 + 1 ,-U J + 0 , 0 1 1 2 tj * A h 1 1 ) / U . 4 

(_• If J(t If If: He 4- ^ ^ CtiMFOlt. ZOAK 

3 2, = />F/USrH 

U (Z. b l'.O.OUl J Gl) I'U 10 
J=i 

ZOrO.UOl 
0 ZO=ZP+0.0149U 

vJ=J+l 

IFCJ.Gt.lZ) GO llJ / 

XF(Z.Gr.ZP) GO Tu b 

1 A = .i 

JZM=J-1 

(;Ht444 4)f 4444 4 CJMl'UTt .iXyc.i., wyotii, aXF, '.'vlF, wX, WK 

1MIjF/V=1 , -alFA 

P F T A = ( Z - 0 . 0 0 1 ) / U . 0 1 4 y - ( 1 A 'I - I ) 

‘iBF i A = 1 . ” B FT A 
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wXOebsOALl-' A^<\il)bLC iwUM , 1Z.W , 1 ) + ht T A*UAbt 
S H 1) K b C I J M , I 4 , i ) + A L. i*' A ♦ u d H, i‘ A ♦ w b cl L< ( I M I J , i Z M , I J + A li h A * b K 1' A ♦ 

S 'i ' ) f<: u ( 1 '<1 J , i Z , 1 ) 

.'I Y uMj = 0 ALl?' A^Urtb'l A* vnu^b C l"'iU!''i / 1 'Zm , Z ) + ob 1 A*OALp' A^ 

$//UKfj ( IMJ VI , lii , ^ ) + AbP’A^L)db l'A*».Uc.bC INU , I'ZM , 2 ) + Abf A^dEl A» 
Sm)Eii( IMj, lZf2) 
wX= J-YP.e'-AiXDtb + UUbE 
« Y = '‘ Yi,)Eb-VKti.e 

UA,X = '-)AuFA*i>l)C,l'A=«‘l^.''i (1''1UH , IZ't, 1 ;_ + bbrA»bAliP A^^UW (IhU'Vs f » 1 ) 

5 + Aijp A*bdETA’*-UW ( J.i‘ U ,lZ.“i , i ) + A bFA * ub X A *Uv'i ( IMLI , 1 i , 1 ) 

Um y = I) a b F a ♦ D d b 1 a b V- ( i M U M , I z Vi , 2 ) + b ET A ♦ U A b F A =*= D w C 1 iM U M , 1 Z , 2 ) 

$ + A j , p A. ♦ D b E r A * i) rt C i Ml , A Zrt , 2 J f A liE A ♦bE X A ♦()« ( A » U , 1 Z , 2 ) 
wX=dX»USrK 
/, Y = %■ y li s I fj 

I ) ‘.' A, = « X F . 

u<*'y = uAy<=F 
bn X t i 

/ WK l i.'E(6, 7 00) 
bu rj i 
Id wy=o. 

AX = US1K*( AbUG(«U6X'H»Z+i .) + 11.2b»AmU + Z)/u,4 
t - ) 3.M I 

dwKllt.t6,7lO) 

/lU I- JK »lAf(5X, • AdU lb UUX iJc nEGlUN’) 

/UU P uK -lA J ( 2X , ' Z IS bAKGc,F. IHA'M O.lS*) 

1 Z = Z + 'JSrK/F 

K {• ■ T M H M 
P, M U 
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APPENDIX 4A 
TABULATED FRONTAL DATA 

Tabulated data for two cold air synoptic fronts are giyen In this 
appendix. The tables have, due to their length, been split Into six parts 
covering columns 1 through 20, 21 through 40, etc. The tabulated values of 
are values relative to the storm motion. The frontal speed ¥ Is given 
at the top of the table for converting W to Its value relative to the 

A ' , . 

ground. 

The case numbers for the two stom fronts studied are listed at the top 
of the table. Also listed at the top of the table are the hor1zonta.l length 
scales for the given wind record. The horizontal extent of each data set 
varies because of the data reduction procedure. Hence, the length of field, 
L, In kilometers and the horizontal grid spacing. Ax, are specified on each 
table. The vertical extent of each field Is taken as 500 m with 50 m 
vertical grid spacing. 
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Table 4A-l(a) Longitudinal Wind Speed, W^, for Cold Front; Case 1, 0307, W^=9o3 m 

L=17„4km, AX=153.7m 
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Tabl e 4A--1 (a ) Conti nued 
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Table 4A-l (b) Lateral Wind Speed, W . for Gold Front; Case 1, 0307, H^=9»3 m 

y ^ 

L=17,4 km, aX= 153„7 m 
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Table 4A-l(b) Continued 
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Table 4A-V(c) Continued 
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Table A4-2(a) Longitudinal Wind Speed, W^, for Cold Front; Case 2, 1050, 
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Table 4A-2(a) Continued 
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Table 4A-2{b) Lateral Wind Speed, W for Cold Front; Case 2, 1050 
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Table 4A-2(c) Vertical Wind Speed* W^, for Cold Front; Case 2, 1050, R =7„1 

X 

L=13.2 km, AX=117.7m 
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APPENDIX 4B 


GRAPHICAL ILLUSTRATION OF FRONTAL WIND SPEED PROFILES 

In this appendix, the wind speed profiles which an aircraft would 
encounter approaching along a 3® glide slope through a cold synoptic front 
are illustrated. Four flight paths having 3® glide slopes and spaced at 
equal increments along the horizontal extent of the data set are shown in 
Figure 4B-1. Figures 4B-2 and 4B-3 illustrate the wind speed profiles along 
the four flight paths. One observes that the winds can be significantly 
different for aircraft approaching at different times through the same front. 
Due to the limited spatial extent of the data set the latter flight paths, 3 
and 4, do not extend to the full 500 m level and consequently are only 
plotted to the height at which th^y pass out of the data set at the right- 
hand margin. 

The streamlines on which the flight paths are drawn are based on wind 
speeds relative to the motion Of the front which is 9.3 m s"^ and 7.1 m s"^ 
for Figures 4B-2 and 4B-3, respectively. The vertical dashed line on the 
longitudinal velocity pi^ofiles illustrate the mean motion of the front, W . 



Case 
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Figure 4B-1 Flight Paths through Streamline Pattern Moving with the Front 
Corresponding to Wind Speed Profiles in Figures 4 B-2 and 4B-3 


-15 -10 -5 0 5 10 15 

( m s ^ ) 


Figure 4B-2 Wind Speed Profiles Along Four Different Flight Paths Through 
Cold Front,, Case 1 











